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EXECUTIVE  SUMMARY 


The  Federal  Aviation  Administration  Technical  Center  owns  and  operates  an  instrumented  B- 
727  aircraft  for  in-house  test  and  assessment  of  runway  friction  and  to  examine  landing  gear 
impact  on  runway  surfaces.  Numerous  strain  gages  were  installed  on  the  landing  gear.  Prior 
analysis  of  the  data  was  concerned  primarily  with  supporting  the  FAA’s  airports  research 
program;  however,  it  has  been  recognized  that  this  data  set  also  has  application  in  correlating 
the  aircraft  gear  loads  to  measured  ground-flight  loads  for  specific  takeoff  and  landing 
conditions  and  runway  operations.  This  in-house  research  affords  the  National  Aging  Aircraft 
Research  Program  (NAARP)  an  opportunity  to  assess  new  data,  independently  obtained, 
which  is  both  of  high  quality  and  unusually  high  sampling  rates  to  examine  the  landing  gear 
loads  for  a  limited  number  of  common  aircraft  maneuvers.  Sixty  sample-per-second  time 
history  traces  were  available  for  each  of  72  events  with  the  following  breakdown:(l)  Takeoff 
analysis,  22  events;  (2)  Landing  analysis,  18  events;  (3)  Runway  exit  analysis,  9  events;  (4) 
Braking  analysis,  6  events;  (5)  S-tum  analysis,  9  events;  and  (6)  Minimum-radius  turn 
analysis,  8  events. 

Multiple  time  trace  data  plots  of  all  maneuvers  illustrate  variable  relationships  as  expected. 
Pitch  rate  was  calculated  for  all  takeoffs  that  had  available  angle  data.  The  average  pitch  rate 
was  calculated  to  be  3.43  degrees  per  second,  while  the  maximum  was  4.63  degrees  per 
second.  Histograms  were  plotted  and  mean,  median,  and  standard  deviations  were  calculated 
for  all  takeoff  and  landing  critical  parameters  although  the  small  number  of  samples  caused  a 
large  amount  of  uncertainty.  Measured  values  were  compared  with  calculated  ones  and  had 
good  agreement,  although  a  great  deal  of  scatter  was  illustrated  for  the  measured  and 
calculated  braking  drag  shear.  This  is  partially  due  to  the  large  bias  in  the  measured  braking 
accelerations  compared  to  the  calculated  accelerations.  Measured  accelerations  agree  well 
with  measured  forces.  Relationships  between  braking  acceleration  and  drag  shear,  lateral 
acceleration  and  side  shear,  and  normal  acceleration  and  vertical  shear  were  determined  for 
the  braking,  exiting,  and  landing  tests  respectively.  A  linear  relationship  exists  between  side 
shear  and  axle  differential  load. 


IX 


L  INTRODUCTION. 


1.1  Background. 

The  uses  and  applications  of  many  of  the  older  commercial  aircraft  operating  today  (i.e., 

DC9,  B-727)  have  slowly  evolved  over  the  past  30  years  since  the  conception  of  the  design. 
Many  common  ground  maneuvers  such  as  aircraft  braking  and  exiting  are  different  due  to 
changes  in  airport  capacity  or  aircraft  use.  Most  of  the  data  available  describing  aircraft  loads 
during  takeoff,  landing,  and  ground  maneuvers  is  as  much  as  30  years  old,  where  prior 
generation  equipment  was  used  acquiring  data  at  possibly  insufficient  data  rates.  This  study 
examines  some  recently  obtained  high  quality,  independent  data  with  high  sample  rates,  to 
measure  and  correlate  landing  gear  internal  and  external  loads  experienced  during  takeoff, 
landing,  and  some  selected  abrupt  ground  maneuvers. 

The  data  will  provide  an  understanding  of  the  relationships  between  the  external  conditions 
(i.e.,  runway  exit  speed)  and  the  internal  loads  (e.g.,  individual  main  gear  loads)  which  may 
be  of  value  in  understanding  service  problems  on  the  older  airframes  and  will  assist  in  the 
development  of  the  design  loads  and  operating  loads  on  new  designs. 


1.2  Scope. 

Data  reduction  of  extensive  ground  and  flight  tests  in  support  of  the  FAA  airport  runway 
friction  research  was  reconfigured  to  allow  for  the  analysis  of  gear  loads.  All  data  were 
acquired  from  the  flight  test  data  acquisition  system  on  the  FAA  Technical  Center’s  research 
transport  B-727  N40  and  most  data  came  from  experiments  performed  to  examine  ground 
handling  performance.  Figure  1  gives  a  plan  view  of  a  B-727- 100QC.  The  flight  and  ground 
tests  in  the  subject  study  include  several  different  maneuvers  such  as  a  takeoff  or  runway  exit. 
Other  data  files  that  have  been  made  available  are  from  tests  to  evaluate  runway  friction,  soft 
arrester  system  foam  performance,  and  other  runway  parameters.  Figure  2  is  a  photograph  of 
the  FAA  Technical  Center’s  B-727  N40  during  taxi.  Appendix  A  provides  tables  of  all  useful 
test  data  files  available  including  parameters  describing  the  test  objectives  and  conditions. 
These  parameters  (i.e.,  gross  weight,  wind  speed,  and  direction)  are  used  for  correlation  of 
measured  internal  loads  for  individual  flight  tests. 

The  purpose  of  this  study  was  to  utilize  existing  high  speed  B-727  data  to  compile  the  flight 
acceleration,  aircraft  system,  and  force  data  for  each  maneuver  into  data  files,  and  to  analyze 
the  flight  information  to  provide  a  typical  picture  of  each  aircraft  maneuver.  Relationships 
describing  takeoff,  landing,  exit,  braking,  S-tum,  and  minimum-radius  turn  aircraft  maneuvers 
were  examined.  These  relationships  provide  some  directly  measured  internal  loads 
information  not  available  in  previous  research  efforts. 
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2.  DATA  ACQUISITION  AND  REDUCTION. 


2.1  Data  Channels. 

Since  the  focus  of  the  testing  was  ground  handling,  most  data  channels  consisted  of  landing 
gear  information  such  as  axle/strut  forces,  wheel  braking  forces,  wheel  steering  angles,  and 
wheel  speeds.  Other  analog  inputs  consist  of  the  standard  control  inputs,  and  engine  pressure 
ratios  (EPRs),  as  well  as  cockpit  and  eg  accelerations.  Most  digital  inputs  came  from  the  air 
data  computer  (ADC)  which  contains  the  inertial  navigation  system’s  (INS)  output  as  well  as 
other  digital  flight  parameters.  These  values  are  typically  flight  path  angles,  headings,  and 
other  data  intrinsic  to  the  aircraft’s  flight  instrument  presentation.  Ground  speed  and  distance 
were  provided  by  a  digital  distance  sensor  which,  by  means  of  a  focused  light  beam,  provided 
a  digital  pulse  every  2.5  mm  of  ground  distance.  A  view  of  this  device  mounted  on  the  front 
landing  gear  can  he  seen  as  figure  3. 

Strain  gages  were  installed  on  all  of  the  landing  gear  axle  assemblies.  The  main  gear 
installation  is  illustrated  in  figure  4.  This  installation  combination  allows  for  the 
measurement  of  vertical  shear,  braking  shear,  and  side  shear.  Although  each  wheel  axle  was 
instrumented  separately,  each  landing  gear  assembly  (left  and  right  wheel  axle)  were 
calibrated  together  by  attaching  the  entire  landing  gear  to  a  load  cell  platform  and 
incrementally  lowering  the  aircraft  onto  the  platform  (see  figure  5).  Although  this  method 
gives  an  accurate  calibration  of  the  three  basic  force  directions,  it  does  not  provide  the 
individual  axle  force  calibration  or  the  cross-coupling  of  the  strain  gage  outputs  from  axle  to 
axle  and  gage  to  gage.  Thus,  the  accuracy  of  axisymmetrical  loads  on  the  landing  gear  has 
not  been  determined.  Appendix  B  describes  the  basic  calibration  procedure  for  the  three 
landing  gears. 


2.2  Data  Acquisition  System  and  Data  Reduction. 

The  data  acquisition  system  installed  on  the  B-727  is  a  Metraplex  hybrid  data  acquisition 
system  which  has  a  series  of  digital  cards  to  acquire  analog  and  different  forms  of  digital  data 
and  creates  a  pulse-code  modulated  (PCM)  encoded  bit  stream  as  output  (see  figure  6).  The 
B-727  data  acquisition  system  has  an  Eidel  PCM  760  decoding  card  to  decode  the  PCM  bit 
stream  from  the  Metraplex  system.  The  PCM  decoding  card  utilizes  two  source  files  to 
decode  the  PCM  bit  stream  and  create  a  data  "frame"  which  contains  all  data  channels,  time 
code,  and  other  important  frame  information.  The  data  frame  is  then  saved  as  a  binary 
sequential  data  file,  along  with  an  ASCII  "playback"  file  which  contains  the  time  code 
information.  The  signal  file  and  the  format  file  are  the  two  source  files  utilized  by  the  EE315 
software  to  make  binary  sequential  data  files  from  the  encoded  PCM  data  stream.  It  is 
important  to  note  that  not  every  test  employed  the  same  signal  file,  or  even  the  same 
channels,  since  the  EE315  software  allows  the  option  of  saving  some  or  all  of  the  channels 
specified  in  the  Metraplex  data  frame. 
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In  order  to  obtain  an  ASCII  data  file  with  a  valid  time  code  and  data  columns  from  the  Eidel 
binary  sequential  file  system,  a  software  package  was  developed  to  allow  the  user  to  convert 
specified  data  from  the  file.  This  software  is  also  used  to  apply  calibration  factors  and 
formulas,  add  or  average  columns,  and  filter  the  data  by  frequency.  This  program  was  written 
in  QBasic  and  allows  the  user  to  obtain  any  data  channels  desired  from  the  data  frame  for  any 
specified  time  frame,  as  well  as  vary  the  output  data  rate.  Because  many  data  channels  were 
redundant  and  unnecessary  for  this  study,  not  all  data  channels  were  converted  to  ASCII  for 
analysis.  Appendix  C  gives  a  list  of  all  data  channels  converted  to  ASCII  with  units  and 
range.  A  listing  of  this  program  can  be  seen  in  appendix  D. 


2.3  Data  Filtering. 

As  is  customary  with  most  analog  force  related  instrumentation  measurements,  the  data  was 
filtered  in  two  stages.  All  strain  gage  signals  are  first  passed  through  a  signal  conditioner, 
which  amplifies  and  then  low-pass  filters  the  data  for  antialiasing.  Most  other  analog  signals 
contain  some  internal  amplification  and  filtering.  The  analog  and  digital  signals  are  then  sent 
to  the  data  acquisition  system  and  stored.  When  the  binary  sequential  data  is  processed  and 
made  into  ASCII  files,  a  fourth-order  Butterworth  HR  digital  filter  is  applied  to  reduce  high 
frequency  noise,  and  to  antialias  signals  that  are  stored  at  a  lower  sample  rate  in  the  ASCII 
files  than  in  the  original  binary  files.  The  effect  of  the  filter  on  ten  seconds  of  data  can  be 
seen  in  figures  7  and  8. 
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3.  ANALYSIS. 


3.1  General  Analysis  Techniques. 

To  reduce  the  data,  the  raw  60  Hz  data  were  digitally  filtered  at  5  Hz  and  decimated  by  2 
(every  second  point  eliminated).  The  cutoff  frequency  (Fc)  was  chosen  to  allow  for  the 
lowest  possible  cutoff  frequency  that  is  significantly  greater  then  the  aircraft  body  natural 
frequency  (3.5  Hz)  and  the  data  were  decimated  to  allow  for  easier  file  management.  This 
gives  a  reduced  data  rate  of  30  Hz  filtered  at  5  Hz.  The  data  were  analyzed  both  in  time 
history  and  for  specific  event  values  such  as  peak  aircraft  longitudinal  acceleration  which  take 
place  at  or  over  some  period  of  time. 

For  analysis  in  time,  certain  time  trace  data  were  correlated  to  verify  the  validity  of  the  data 
and  also  to  illustrate  certain  key  parameter  relationships.  Correlation  in  time  of  two 
parameters  was  achieved  by  plotting  each  parameter  on  a  separate  axis  for  a  selected  time 
period.  Relationships  between  the  two  parameters  could  then  be  observed  and  analyzed. 

When  relating  two  parameters  in  time,  they  were  related  over  a  specific  time  period  pertinent 
to  the  maneuver,  and  only  data  obtained  during  this  specific  time  period  were  considered  for 
correlation.  All  leading  and  lagging  data  points  were  "trimmed"  to  reduce  data  scatter.  These 
special  trimmed  files  also  provided  a  basis  for  determining  specific  event  values,  which  were 
often  an  average  value  over  a  small  interval  of  time  during  the  event. 

For  the  analysis  of  specific  event  parameters,  each  aircraft  maneuver  specific  event  value  was 
determined  and  tabulated.  Each  specific  event  parameter  was  chosen  for  its  relevance  to  each 
maneuver  (i.e.,  landing  normal  acceleration  at  touchdown).  These  values  were  tabulated  for 
evaluation  and  possible  correlation  between  tests,  whereas  each  flight  segment  had  a  prime 
focus  for  analysis.  Statistical  analysis  of  these  critical  values  was  also  performed  for  the 
takeoff  and  landing  events,  where  a  histogram  of  each  specific  event  parameter  was  plotted 
and  mean,  median,  and  standard  deviations  were  calculated  for  each  histogram  using  CoStat,  a 
commercial  software  package.  Specific  event  parameters  were  also  compared  with  calculated 
values  to  gage  the  accuracy  of  the  data  and  to  verify  the  means  by  which  acceleration  is 
translated  to  force.  Specific  event  values  were  also  correlated  with  each  other  to  determine 
possible  relationships,  although  correlation  of  these  values  was  difficult  due  to  the  limited 
amount  of  test  events  for  each  maneuver.  However,  even  with  this  limited  data,  some 
significant  relationships  were  determined. 


3.2  Takeoff  Data  Analysis. 

The  primary  focus  of  the  takeoff  data  analysis  was  on  takeoff  velocity  and  rotation  data. 

Some  takeoff  tests  did  not  have  aircraft  flight  path  angles  available  for  analysis,  thus  reducing 
the  number  of  flight  events  available  for  analysis.  The  key  parameters  tabulated  for  each 
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takeoff  test  other  than  the  basic  test  information  were  peak  longitudinal  acceleration  (Nx), 
maximum  ground  speed  (Vmax),  runway  distance  used  (AD),  time  to  liftoff  (tlo),  maximum 
pitch  angle  (0max),  and  maximum  pitch  rate  (PR^).  The  maximum  takeoff  velocity  adjusted 
for  wind  speed  was  also  tabulated.  These  parameters  are  illustrated  in  figure  9.  The  peak 
longitudinal  acceleration  was  approximated  by  drawing  a  best  fit  curve  to  the  acceleration 
time  trace,  and  then  determining  the  peak  value.  The  ground  speed  and  distance  sensor  was 
located  on  the  nose  gear  (see  section  2.2)  and  these  values  had  to  be  estimated  after  nose  gear 
liftoff.  Maximum  ground  speed  was  assumed  to  be  the  ground  speed  at  the  start  of  rotation, 
and  the  runway  distance  utilized  after  rotation  started  was  approximated  by  dividing  the 
maximum  ground  speed  by  the  time  remaining  until  main  gear  liftoff.  Histograms  of  each 
critical  parameter  were  calculated  and  plotted  for  all  tests  and  the  mean  value,  median  value, 
and  standard  deviation  were  calculated.  Also,  takeoff  speed  was  correlated  with  takeoff 
distance  and  time  to  liftoff  to  determine  the  relationship  between  these  takeoff  parameters. 


3.3  Landing  Data  Analysis. 

The  primary  focus  of  the  landing  data  analysis  was  data  available  on  vertical  eg  acceleration 
and  peak  axle  vertical  shear  at  aircraft  touchdown.  The  key  parameters  tabulated  for  each 
landing  other  than  the  basic  flight  information  were  peak  vertical  acceleration  (Nz),  nose  gear 
touch  down  ground  speed  (Vtd),  touchdown  pitch  angle  (0td),  and  peak  main  landing  gear 
vertical  shear  (VSMLG).  The  touchdown  ground  speed  adjusted  for  wind  velocity  was  also 
tabulated.  These  parameters  are  illustrated  in  figure  10.  Histograms  were  calculated  and 
plotted  for  each  critical  parameter  and  the  mean  value,  median  value,  and  standard  deviation 
were  calculated.  In  an  attempt  to  gage  the  accuracy  of  the  data,  landing  vertical  shear  was 
calculated  using  the  following  equation  and  compared  to  the  measured  vertical  shear. 

"W  =  Nt  WAC  (1) 


Where  WAC  is  the  aircraft  landing  weight,  vertical  shear  and  normal  acceleration  are  single 
peak  values.  Peak  landing  gear  vertical  shear  was  correlated  with  normal  acceleration  and 
also  pitch  angle  at  touchdown  was  correlated  with  aircraft  weight  in  an  attempt  to  determine 
the  relationship  between  these  landing  parameters. 


3.4  Runway  Exit  Data  Analysis. 

The  main  focus  of  the  exit  analysis  was  to  obtain  the  relationship  between  landing  gear  side 
force  and  lateral  acceleration,  as  it  relates  to  different  exit  speeds  on  a  standard  runway  exit 
ramp.  Measured  side  shear  was  correlated  in  time  with  measured  lateral  acceleration  for  a 
40-,  50-,  and  60-knot  exit  to  illustrate  the  relationship  between  the  two  parameters  and  verify 
the  quality  of  the  measurements.  As  stated  in  section  3.1,  only  the  period  of  time  during 
which  the  aircraft  was  actually  exiting  the  runway  was  analyzed.  Time  averages  of  some 
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critical  values  were  utilized  to  allow  for  correlation  of  tabulated  critical  parameters.  Key 
parameters  tabulated  were  exit  time  interval  (At),  exit  distance  interval  (AD),  average  ground 
speed  (Vave),  average  lateral  acceleration  (Ny),  average  main  landing  gear  side  shear  (SSMLG), 
and  peak  main  landing  gear  vertical  shear  (VSMLG).  These  values  are  illustrated  in  figure  11. 
To  verify  measured  average  lateral  acceleration  and  average  side  shear,  values  were  calculated 
utilizing  other  measured  values.  To  calculate  average  acceleration,  the  average  turn  radius 
was  first  calculated. 

-  A  D  (2) 

r  =  - 

0 


In  the  above  equation  0  is  the  angle  representing  the  change  in  path  of  the  aircraft  and  the 
over  bar  represents  a  time  averaged  value.  Assuming  the  lateral  acceleration  is  equal  to  the 
centrifugal  radial  force,  the  following  equation  was  used  to  calculated  average  lateral 
acceleration. 


(3) 


Again,  the  over  bar  represents  a  time  averaged  value.  Average  side  shear  can  be  calculated 
using  average  measured  lateral  acceleration  (Ny)  and  aircraft  weight  (WAC). 


SSMLG  =  Xy  WM 


(4) 


These  simple  calculations  can  help  validate  the  accuracy  of  some  measured  parameters.  To 
compare  side  force  with  lateral  acceleration  in  time,  side  force  coefficient  was  calculated 
using  the  following  formula. 


S^tot 

“  W 

”  AC 


(5) 


In  the  above  equation  Css  is  a  local  time  value  at  the  local  total  side  shear  SSTOT.  Note  that 
the  above  equation  represents  the  reverse  calculation  of  equation  (4).  Also,  lateral 
acceleration  was  correlated  with  exit  speed  and  landing  gear  side  force.  Side  Shear  was  also 
correlated  with  exit  speed  to  illustrate  this  relationship.  The  aircraft  wheel  brakes  were  not 
applied  during  the  exit  tests  while  the  aircraft  was  traveling  through  the  turn.  All  the  data 
presented  for  the  exit  tests  were  measured  with  the  wheels  freely  rolling.  The  optional  nose 
wheel  brakes  were  disabled. 


3.5  Braking  Data  Analysis. 


Braking  data  were  analyzed  in  an  attempt  to  obtain  a  relationship  between  braking 
acceleration  and  axle  drag  shear.  Measured  drag  shear  was  correlated  in  time  with  measured 
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longitudinal  acceleration  for  several  braking  events  of  different  intensities  to  illustrate  the 
relationship  between  the  two  parameters  and  to  verify  the  quality  of  the  measurements.  Only 
the  portion  of  each  braking  test  which  involved  actual  braking  was  analyzed.  The  key 
parameters  utilized  for  each  of  these  braking  events  included  change  in  speed  (AV),  time  to 
change  speed  (At),  braking  distance  (AD),  average  braking  acceleration  (Nx),  average  total 
drag  shear  (DSMLG),  and  maximum  braking  pressure  (Pbrk).  Also,  the  initial  ground  speed  at 
the  start  of  braking  was  recorded  (Vstrt).  These  values  are  illustrated  in  figure  12.  To  verify 
measured  average  acceleration  and  average  drag  shear,  values  were  calculated  utilizing  other 
measured  values.  Average  acceleration  can  be  calculated  by  using  change  in  velocity  (AV) 
and  time  (At)  using  the  following  formula. 


AV 

At 


(6) 


Average  drag  shear  can  be  calculated  using  average  measured  acceleration  (Nx)  and  aircraft 
weight  (WAC). 


DSmlg  =  Nx  WAC 


(7) 


The  average  acceleration  tabulated  was  correlated  with  the  average  drag  force  tabulated  to 
obtain  a  relationship  between  the  two  parameters.  All  braking  tests  utilized  were  performed 
on  dry  pavement  except  TST71001.  Note,  as  for  the  exit  tests,  the  nose  wheel  brakes  were 
disabled. 


3.6  S-Turn  Data  Analysis. 

The  primary  focus  of  S-turn  data  analysis  was  to  examine  the  vertical  shear  axle  differential 
loads  experienced  during  lateral  acceleration  due  to  the  aircraft’s  lateral  movement.  Axle 
differential  load  is  defined  as  the  difference  between  the  left  axle  vertical  shear  and  the  right 
axle  vertical  shear  for  a  left  turn.  Measured  side  shear  was  correlated  in  time  with  measured 
lateral  acceleration  for  a  40-,  60-,  and  80-knot  S-turn  to  illustrate  the  relationship  between  the 
two  parameters  and  verify  the  quality  of  the  measurements.  Also  to  verify  the  validity  of  the 
data,  the  side  force  coefficient  was  plotted  with  lateral  acceleration  in  time.  As  stated  in 
section  3.1,  only  the  period  of  time  during  which  the  aircraft  was  actually  executing  an  S-tum 
was  analyzed.  Time  averages  of  some  critical  values  were  utilized  to  allow  for  correlation  of 
tabulated  critical  parameters.  Key  parameters  tabulated  were  S-tum  start  speed  (Vstrt),  time 
interval  (At),  average  lateral  acceleration  (Ny),  average  landing  gear  side  shear  for  each  main 
landing  gear  (SSLMG,  SSRMG),  and  average  landing  gear  vertical  axle  differential  load 
(ADLrmg,  ADLlMG)  for  each  main  landing  gear.  The  maximum  main  landing  gear  vertical 
axle  differential  load  (ADLmlg)  was  also  tabulated.  These  parameters  are  illustrated  in  figure 
13.  To  verify  measured  average  side  shear,  values  were  calculated  using  measured  average 
lateral  acceleration  and  aircraft  weight. 
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(8) 


SSMLG  =  "y  WAC 


In  this  equation  SSMLG  is  the  total  side  shear  on  both  landing  gear.  Side  shear  was  correlated 
with  axle  differential  load  and  also  with  lateral  acceleration  in  an  attempt  to  determine  the 
relationship  between  these  parameters.  Braking  was  not  used  during  the  S-Turn  tests  and  the 
optional  nose  wheel  brakes  were  disabled  for  the  testing. 

3.7  Minimum-Radius  Turn  Data  Analysis. 

The  primary  focus  of  minimum-radius  turn  data  analysis  was  the  examination  of  the  nose  gear 
side  shear  and  vertical  shear  axle  differential  loads  experienced  during  the  maneuver. 

Measured  side  shear  on  the  nose  gear  was  plotted  in  time  with  measured  nose  gear  axle 
differential  load  for  a  clockwise  (CW)  and  counter  clock  wise  (CCW)  minimum-radius  turn  to 
illustrate  the  relationship  between  the  two  parameters.  An  arbitrary  ten-second  period  of  time 
at  which  the  aircraft  was  actually  executing  the  turn  was  analyzed  when  considering  time 
correlation  and  averages.  Key  parameters  tabulated  were  minimum-radius  turn  time  interval 
(At),  average  nose  gear  tangential  velocity  (Vtan),  average  nose  gear  side  shear  (SSNG),  and 
average  nose  gear  vertical  axle  differential  load  (ADL^).  The  maximum  nose  gear  vertical 
axle  differential  load  (ADL^)  was  also  tabulated.  These  parameters  are  illustrated  in  figure 
14.  Nose  gear  side  shear  was  correlated  with  axle  differential  load  in  an  attempt  to  determine 
the  relationship  between  these  parameters.  Braking  was  not  used  during  the  minimum  radius 
turn  tests  and  the  optional  nose  wheel  brakes  were  disabled  for  the  testing. 
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4.  RESULTS. 


4.1  Takeoffs. 

Figures  15  and  16  illustrate  aircraft  acceleration,  ground  speed,  and  aircraft  angle  data  for  two 
typical  takeoffs  of  the  instrumented  B-727  N40.  Note  the  peak  longitudinal  acceleration 
occurred  early  during  the  takeoff  run  time  trace.  The  start  of  rotation  is  illustrated  by  the 
ground  speed  going  to  zero  at  approximately  the  same  time  the  pitch  angle  reaches  zero.  Roll 
angle  changes  very  little  during  the  takeoff  run.  Liftoff  is  signified  by  a  drastic  decrease  in 
noise  on  the  normal  acceleration  data  trace.  Pitch  data  for  each  of  the  illustrated  aircraft  tests 
are  given  in  figures  17a  and  17b.  Maximum  pitch  rate  during  TST63006  (2.57  deg/s)  was 
much  lower  than  in  TST11408  (4.17  deg/s). 

Table  1  contains  the  critical  takeoff  parameters  discussed  in  section  3.2.  The  cells  marked 
with  an  "X"  indicate  information  that  was  not  available  for  one  or  more  reasons.  The  first 
four  columns  are  information  relative  to  the  flight  test  itself  (aircraft  weight,  wind  speed, 
etc.),  while  the  other  seven  columns  represent  the  critical  parameters. 

To  determine  if  the  critical  parameters  collected  are  a  good  measure  of  expected  values  for  an 
aircraft  of  similar  or  same  type,  a  histogram  of  each  critical  parameter  was  plotted  to 
determine  if  the  scatter  in  the  data  are  indeed  a  random  deviation  from  an  average  value. 
Figures  18  through  24  illustrate  histograms  of  the  critical  takeoff  parameters,  with  each 
histogram  displaying  the  sample  size,  mean,  median,  and  standard  deviation  of  the  data  field. 
Most  histograms  have  a  median  value  close  to  the  mean  value,  which  illustrates  good 
symmetry. 

It  would  be  desirable  if  the  takeoff  distance  could  be  related  to  takeoff  velocity  and  aircraft 
weight.  Figure  25  gives  a  correlation  plot  of  takeoff  distance  versus  takeoff  speed  to 
illustrate  a  potential  relationship.  Although  some  data  scatter  exists,  a  linear  relationship  can 
be  observed.  Figure  26  gives  a  correlation  plot  of  takeoff  distance  versus  aircraft  weight. 

The  plot  illustrates  a  large  amount  of  data  scatter. 

4.2  Landings. 

Figure  27  and  28  illustrate  aircraft  accelerations,  landing  gear  vertical  shear,  and  pitch  angle 
data  for  two  typical  landings  of  the  instrumented  B-727  N40.  Note  the  maximum  normal 
acceleration  did  not  always  occur  at  touchdown,  although  it  does  for  the  two  given  landings 
and  it  is  this  touchdown  normal  acceleration  which  is  the  focus  of  the  landing  acceleration 
data.  Touchdown  has  occurred  when  there  exists  a  sudden  increase  in  the  vertical  shear  and 
the  drastic  increase  in  noise  on  the  normal  and  longitudinal  acceleration  data  trace.  The 
landing  pitch  angle  time  trace  is  different  for  each  individual  landing,  however  the  maximum 
pitch  angle  at  touchdown  is  a  fairly  consistent  value. 
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Table  2  contains  the  critical  landing  parameters  discussed  in  section  3.3.  The  cells  marked 
with  an  "X"  indicate  information  that  was  not  available  for  one  or  more  of  the  reasons 
discussed  m  section  2.3.  The  first  four  columns  are  information  relative  to  the  flight  test 

itself  (aircraft  weight,  wind  speed,  etc.),  while  the  other  five  columns  represent  the  measured 
landing  parameters. 

As  discussed  in  section  3.2,  the  peak  vertical  shear  values  measured  at  touchdown  were 
compared  to  values  calculated  with  acceleration  and  aircraft  weight.  This  comparison  is 
shown  graphically  as  figure  29.  Although  agreement  is  not  always  good,  many  calculated 
values  have  less  than  10  percent  error  from  the  measured.  Acceleration  and  weight 
consequently  are  not  the  only  two  parameters  which  correlate  with  shear.  There  is  not 
sufficient  data  in  this  study  to  perform  a  more  in  depth  analysis,  however  given  that  peak 
vert, cal  shear  is  a  function  of  N,  and  W«  (Equation  4),  vertical  shear  is  predicted  with  a 
correlation  coefficient  of  0.798.  This  is  likely  explained  due  to  ground  friction  forces  and 
airplane  dynamics  affecting  the  aircraft  eg  accelerations. 

To  determine  if  the  critical  parameters  collected  are  a  good  measure  of  expected  values  for  an 
aircraft  of  similar  or  same  type,  a  histogram  of  each  critical  parameter  was  plotted  to 
determine  if  the  scatter  in  the  data  is  indeed  random  deviation  from  an  average  value. 

Figures  30  through  34  illustrate  histograms  of  the  critical  landing  parameters,  with  each 
histogram  displaying  the  data  fields  size,  mean,  median,  and  standard  deviation.  Most  median 
values  are  close  to  the  mean,  indicating  symmetry. 

It  would  be  desirable  to  correlate  touchdown  vertical  shear  with  normal  acceleration 
experienced  Figure  35  gives  a  correlation  plot  of  touchdown  vertical  shear  versus  touchdown 
normal  acceleration  with  a  correlation  coefficient  of  0.784.  Although  a  relationship  is 
evident,  the  problems  with  relating  these  two  parameters  discussed  in  the  previous  paragraph 
are  illustrated  by  the  data  scatter.  Figure  36  gives  a  correlation  plot  of  touchdown  pitch  angle 
versus  aircraft  weight  with  a  correlation  coefficient  of  0.512.  Again,  a  relationship  is  evident 

correlation^  Confldence  is  illustrated  by  the  large  amount  of  data  scatter  causing  a  poor 


4.3  Runway  Exits. 


igures  37,  38,  and  39  illustrate  aircraft  acceleration,  landing  gear  vertical  shear,  and  landing 
gear  side  shear  for  three  separate  runway  exits  at  40,  50,  and  60  knots  respectively.  Viewing 
the  succession  of  graphs  illustrates  the  increase  in  lateral  acceleration  and  axle  shears  as  the 
exit  speed  increases.  To  determine  if  side  shear  magnitude  and  time  behavior  are  accurate 
and  reasonable,  the  side  force  coefficient  was  calculated  and  plotted  along  side  lateral 
acceleration.  Figure  40  gives  the  side  force  coefficient  comparison  which  illustrates  a 
tendency  for  the  side  force  coefficient  to  be  high,  particularly  at  heavy  side  loading  This 
could  be  due  to  torque  on  the  strut  amplifying  the  side  force,  or  to  absorption  of  force  by  the 
aircraft  suspension,  giving  less  acceleration  at  the  aircraft  eg.  A  biased  high  calibration  factor 
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would  also  cause  this  phenomena.  Landing  gear  side  shear  was  plotted  against  lateral 
acceleration  in  time  for  the  three  exits  illustrated,  and  presented  as  figure  41.  As  mentioned 
in  section  3.4  only  the  time  values  during  the  aircraft  exit  were  used  in  making  the  graph. 
Although  the  plot  exhibits  modest  scatter,  a  good  linear  relationship  is  observed.  The  weight 
of  the  aircraft  varied  from  134,726  pounds  to  132,326  pounds,  accounting  for  some  of  the 
plot  scatter,  while  the  least  scatter  is  observed  at  high  acceleration  values  with  the  side  shear 
peaking  at  about  32,000  pounds. 

Table  3  contains  the  critical  exit  parameters  discussed  in  section  3.4.  The  first  two  columns 
are  information  relative  to  the  flight  test  itself,  while  the  other  five  columns  represent  the 
critical  parameters.  As  discussed  in  section  3.4  most  of  these  values  are  time  averages, 
allowing  for  easy  verification  of  the  validity  of  the  data.  To  help  verify  measured  lateral 
accelerations,  a  graph  comparing  measured  and  calculated  lateral  accelerations  is  contained  in 
figure  42.  Note  the  excellent  agreement  in  the  regular  and  high  speed  exits.  A  comparison 
of  measured  and  calculated  side  shear  is  plotted  in  figure  43.  Again  excellent  agreement  can 
be  observed  between  measured  and  calculated  values,  illustrating  the  high  quality  of  the  data. 

Figure  44  is  a  correlation  plot  of  average  exit  speed  versus  average  lateral  acceleration 
tabulated  in  the  critical  value  tables.  The  linear  relationship  observed  illustrates  the 
consistency  of  lateral  accelerations  observed  on  any  given  exit.  To  determine  how  side  shear 
relates  to  lateral  acceleration,  a  correlation  plot  of  the  two  parameters  was  plotted  as  figure 
45.  Similar  good  linearity  is  observed  as  is  expected.  The  lateral  acceleration  critical 
parameter  is  an  average  based  on  the  aircraft  velocity  during  the  exit,  which  is  not  subject  to 
the  ground  roughness  error  seen  in  peak  acceleration  data.  Figure  46  is  a  cross  correlation  of 
the  previous  two  figures  giving  the  relationship  between  average  measured  side  shear  and  exit 
speed. 

4.4  Braking. 

Figures  47,  48,  and  49  illustrate  aircraft  acceleration,  ground  speed,  and  landing  gear  drag 
shear  for  three  separate  braking  tests  with  heavy,  normal,  and  light  braking  respectively.  The 
heavy  braking  test  was  run  with  maximum  braking  effort  applied  by  the  pilot,  and  the  anti¬ 
skid  system  was  in  operation  during  the  test  (  a  Mark  II  Hydro-Aire  antiskid  system  was 
installed  on  the  test  aircraft).  Viewing  the  succession  of  graphs  illustrates  the  decrease  in 
acceleration  and  axle  drag  shears  as  the  braking  intensity  decreases.  To  determine  the  drag 
shear  loading  behavior  of  the  landing  gear,  drag  shear  was  plotted  against  longitudinal 
acceleration  for  the  three  braking  tests  illustrated  in  figure  50.  As  mentioned  in  section  3.5 
only  the  values  during  which  the  aircraft  was  braking  were  used.  The  data  exhibits  a 
significant  scatter,  as  well  as  displaying  different  data  grouping  slopes.  A  possible 
explanation  for  the  different  slope  of  the  heavy  breaking  data  may  be  dynamic  effects  due  to 
the  antiskid  device  used  in  the  B-727  braking  system. 

Table  4  contains  the  critical  braking  parameters  discussed  in  section  3.5.  The  first  two 
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columns  are  information  relative  to  the  flight  test  itself,  while  the  other  seven  columns 
represent  the  critical  parameters.  As  discussed  in  section  3.5  most  of  these  values  are  time 
averages,  allowing  for  easy  verification  of  the  validity  of  the  data.  To  help  verify  measured 
longitudinal  accelerations,  a  graph  comparing  measured  and  calculated  longitudinal 
accelerations  is  contained  in  figure  51.  The  calculated  values  are  consistently  biased  high  as 
compared  to  the  measured  values.  A  comparison  of  measured  and  calculated  drag  shear  is 
plotted  as  figure  52.  This  plot  exhibits  normal  scatter  and  is  not  biased  towards  either 
calculated  or  measured  data. 

Figure  53  is  a  correlation  plot  of  average  longitudinal  acceleration  versus  the  average 
measured  drag  shear.  The  poor  correlation  observed  is  probably  a  reflection  of  the 
inconsistencies  observed  in  figure  50  discussed  in  the  first  paragraph  in  this  section. 


4.5  S-Turns. 

Figures  54,  55,  and  56  illustrate  aircraft  acceleration,  ground  speed,  and  landing  gear  drag 
shear  for  three  separate  S-turns  of  80,  60,  and  40  knots  respectively.  Although  the  slower 
speed  turns  reflect  the  longer  time  needed  to  travel  the  same  distance  at  a  lower  velocity,  the 
acceleration,  vertical  shears,  and  side  shears  appear  about  the  same.  This  is  most  likely  due 
to  pilot  input  based  on  his  understanding  of  the  capability  of  the  aircraft.  To  determine  if 
side  shear  magnitude  and  time  behavior  are  accurate  and  reasonable,  side  force  coefficient 
was  calculated  and  plotted  along  side  lateral  acceleration.  Figure  57  illustrates  the  side  force 
coefficient  comparison  which  illustrated  a  tendency  for  the  side  force  coefficient  to  be  high, 
particularly  at  heavy  side  loading.  This  could  be  due  to  torque  on  the  strut  amplifying  the 
side  force,  or  to  absorption  of  force  by  the  aircraft  suspension,  giving  less  acceleration  at  the 
aircraft  eg.  A  biased  high  calibration  factor  would  also  cause  this  phenomena.  Figures  58, 

59,  and  60  are  graphs  of  side  shear  and  axle  differential  load  plotted  in  time.  With  the 
exception  of  the  left  main  gear  of  TST63008,  side  shear  tends  to  relate  1:1  with  axle 
differential  load.  A  basic  force  analysis  of  a  landing  gear  under  lateral  acceleration  tends  to 
agree  with  this  result.  Appendix  E  gives  a  loading  diagram  and  analysis  which  illustrates  that 
the  relationship  between  side  shear  and  axle  differential  load  is  based  on  the  ratio  of  twice  the 
loaded  tire  radius  to  the  lateral  patch  distance,  which  is  approximately  1.2.  To  determine  the 
side  shear  loading  behavior  of  the  landing  gear,  side  shear  was  plotted  in  time  against  lateral 
acceleration  for  the  three  S-turn  tests  illustrated  in  figure  61.  As  mentioned  in  section  3.1 
only  a  small  portion  of  the  data  obtained  during  the  time  the  aircraft  was  executing  a  turn 
were  used.  The  data  exhibits  an  excellent  linear  relationship.  Also,  side  shear  was  plotted  in 
time  against  axle  differential  load  and  can  be  seen  as  figure  62.  This  graph  also  displays 
little  scatter  and  a  good  linear  relationship,  although  distinct  "runners"  of  data  protrude 
vertically  from  the  data  scatter,  illustrating  that  another  mechanism  may  cause  changes  in  axle 
load  differential  while  side  shear  remains  constant. 

Table  5  contains  the  critical  S-turn  parameters  discussed  in  section  3.6.  The  first  two 
columns  are  information  relative  to  the  flight  test  itself,  while  the  other  six  columns  represent 
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the  critical  parameters.  As  discussed  in  section  3.6  most  of  these  values  are  time  averages, 
allowing  for  easy  verification  of  the  validity  of  the  data.  To  help  verify  measured  side  shear, 
a  graph  comparing  measured  and  calculated  side  shear  is  contained  in  figure  63.  The 
measured  values  are  consistently  biased  high  as  compared  to  the  values  calculated  with 
measured  acceleration.  This  agrees  with  figure  57,  which  also  illustrates  a  high  biased  toward 
the  measured  values. 

Figure  64  is  a  correlation  plot  of  average  lateral  acceleration  versus  the  average  measured  side 
shear.  These  short  time  average  values  illustrate  a  good  linear  relation.  Figure  65  is  a 
correlation  plot  of  average  measured  side  shear  and  average  measured  axle  differential  load. 
Again  these  values  relate  well  linearly. 


4.6  Minimum-Radius  Turns. 

Figures  66  and  67  illustrate  nose  gear  ground  speed,  nose  gear  vertical  shears,  and  nose  gear 
side  shear.  The  trend  of  the  side  shear  traces  is  for  the  side  shear  to  be  approximately  zero 
when  the  aircraft  is  stationary  and  to  increase  to  a  value  greater  than  one  half  the  vertical 
shear  force  as  the  aircraft  begins  to  move.  The  increase  occurs  over  a  short  distance,  when 
speed  is  low  and  remains  fairly  constant  throughout  the  test.  The  difference  between  the  left 
and  right  vertical  shears  shows  the  same  trend,  indicating  that  the  primary  cause  of  the 
differential  vertical  shears  is  rotation  of  the  gear  due  to  the  side  shear. 

Video  tape  of  the  nose  gear  taken  during  the  minimum-radius  turns  clearly  shows,  from  the 
deformation  of  the  tires,  the  presence  of  significant  side  shears  and  significant  differential 
vertical  loads,  confirming  the  validity  of  the  data  in  figures  66  and  67.  Figure  68  is  a  single 
frame  taken  from  the  video  tape  and  shows  the  nose  gear  traveling  toward  the  camera. 

Since  the  magnitude  of  the  side  shears  shows  little  variation  with  speed  once  the  aircraft  is 
moving,  it  appears  that  the  nose  gear  tire  side  forces  are  generated  to  resist  the  moment  about 
the  eg  from  side  forces  and  moments  at  the  main  gear  generated  due  to  scrubbing  and  not  to 
resist  the  centrifugal  force  at  the  eg.  In  this  connection  it  should  be  noted  that  the  speed 
traces  shown  in  figures  66  and  67  are  for  the  speed  of  the  nose  gear  measured  perpendicular 
to  the  nose  gear  axle.  Speed  at  the  eg  was  not  measured,  but  from  the  ratio  of  eg  and  nose 
gear  distances  to  the  center  of  turn,  it  would  be  approximately  one  fourth  the  speed  of  the 
nose  gear.  At  the  maximum  nose  gear  speed  shown  in  the  traces  of  12  mph,  a  maximum 
lateral  acceleration  at  the  eg  of  approximately  1.6  ft/s2  (0.05  g)  is  obtained.  That  there  was 
significant  scrubbing  at  the  main  gear  wheels  on  the  inside  of  the  turn  was  shown  by  dark 
circular  rubber  tracks  deposited  on  the  pavement  by  the  tires  (see  figure  69).  The  radius  of 
the  tire  tracks  showed  that  the  center  of  the  turn  was  about  3.5  feet  outside  the  center  of  the 
main  gear  on  the  inside  of  the  turn,  giving  a  turn  radius  to  the  eg  of  approximately  13  feet. 

Figure  70a  and  70b  illustrate  nose  gear  side  shear  as  compared  to  nose  gear  vertical  axle 
differential  load  in  time.  Although  there  are  discrepancies  in  both  tests  between  the  side 
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shear  and  axle  differential  load,  these  can  be  generally  associated  with  changes  in  speed.  In 
test34  differential  load  and  side  shear  agree  very  well  during  the  first  45  seconds  of  the  test 
when  the  nose  wheel  speed  is  increasing.  Figure  71  is  a  correlation  in  time  of  nose  gear  side 
shear  and  axle  differential  load  for  both  tests  discussed  above  for  an  arbitrary  10-second  time 
interval.  Although  the  data  from  each  test  correspond  very  well,  different  characteristics  are 
shown  for  clockwise  and  counter-clockwise  turns. 

Table  6  contains  the  critical  minimum-radius  turn  parameters  discussed  in  section  3.7.  The 
first  three  columns  are  information  relative  to  the  flight  test  itself,  while  the  other  five 
columns  represent  the  critical  parameters.  As  discussed  in  section  3.7  most  of  these  values 
are  time  averages,  made  over  an  arbitrary  10-second  period  during  the  turn,  allowing  for  easy 
verification  of  the  validity  of  the  data. 

Figure  72  illustrates  a  correlation  of  nose  gear  average  side  shear  and  average  axle  differential 
load.  This  graph  illustrates  two  distinct  groupings  for  the  clockwise  and  the  counter¬ 
clockwise  turns  which  agree  with  figure  7 1 .  The  cause  of  this  behavior  has  not  been 
determined. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS. 


From  this  study  of  flight-ground  loads  measured  on  the  B-727  N40,  the  following  conclusions 
can  be  drawn. 

o  Takeoff  data  had  measured  accelerations  consistent  with  known  reality  and 
calculated  maximum  pitch  rate  had  an  average  of  3.43  deg/s  for  the  examined 
takeoffs  with  the  maximum  pitch  rate  observed  being  4.17  deg/s.  Although 
histograms  were  poor,  mean  values  were  consistent  with  median  values.  Most 
attempts  at  correlating  specific  event  values  illustrated  a  linear  relationship, 
although  correlation  plots  exhibit  a  great  deal  of  scatter. 

o  Landing  data  were  consistent  with  known  reality  and  measured  specific 
event  values  for  vertical  shear  were  consistent  with  measured  vertical 
acceleration  values.  Although  histograms  were  poor,  mean  values  were 
consistent  with  determined  median  values.  Specific  event  vertical  shear  values 
correlated  fairly  well  with  vertical  acceleration  and  aircraft  weight,  but 
touchdown  pitch  angle  correlated  poorly. 

o  Runway  exit  time  trace  data  illustrated  the  expected  force-acceleration 
relationships,  and  measured  side  shear  correlated  very  well  with  measured 
accelerations.  Measured  specific  event  values  were  consistent  with  calculated 
values  and  a  direct  relationship  was  established  between  exit  speed  and 
measured  lateral  accelerations.  Side  Shear  force  correlates  well  with  exit  speed 
and  lateral  acceleration. 

o'  Aircraft  braking  time  trace  data  illustrated  the  expected  force-acceleration 
relationships,  and  measured  drag  shear  correlated  well  with  measured 
longitudinal  accelerations  with  different  relationships  being  observed  for  the 
different  levels  of  braking.  Measured  specific  event  values  were  consistent 
with  calculated  values,  with  the  measured  accelerations  being  consistently 
higher  than  the  calculated.  Attempts  to  correlate  drag  force  with  longitudinal 
acceleration  were  difficult  due  to  a  small  number  of  tests. 

o  S-tum  time  trace  data  were  consistent  with  known  reality,  and  measured 
lateral  acceleration  agreed  well  with  measured  side  shear.  Main  landing  gear 
side  shear  and  axle  differential  load  illustrate  a  distinct  relationship  which 
relates  linearly  in  time.  Specific  event  values  for  lateral  acceleration,  main 
landing  gear  side  shear,  and  axle  differential  load  all  correlate  well,  illustrating 
their  interwoven  relationship. 

o  Minimum-radius  turn  time  trace  data  was  consistent  with  known  reality,  and 
measured  nose  gear  side  shear  related  well  with  measured  nose  gear  axle 
differential  load,  although  it  is  unclear  how  they  relate  in  time.  Specific  event 
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values  for  nose  gear  side  shear  and  axle  differential  load  agree  poorly. 


Additional  takeoff  and  landing  test  events  would  be  useful  in  obtaining  better  histograms,  and 
therefore  more  accurate  statistical  information.  Similarly,  more  runway  exit  data  with  many 
different  exit  velocities  would  be  useful  to  allow  for  more  complete  correlation  plots  between 
landing  gear  side  shear,  lateral  acceleration,  and  exit  velocity.  More  consistent  braking  data 
would  be  useful  in  more  accurately  predicting  drag  shear.  Multiple  braking  runs  with  and 
without  the  antiskid  device  in  use  could  also  help  determine  why  different  braking  slopes  give 
different  slope  scatter  patterns  on  a  correlation  plot  of  longitudinal  acceleration  versus  drag 
force.  Although  an  obvious  relationship  between  side  shear  and  axle  differential  load  exists,  a 
recalibration  of  the  B-727  landing  gear  would  be  needed  to  obtain  the  definitive  relationship. 
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FIGURE  1.  DIAGRAM  OF  B-727  ILLUSTRATING  GEOMETRY  AND  SIGNIFICANT  DISTANCES 


FIGURE  5.  PHOTO  OF  LANDING  GEAR  CALIBRATION  ILLUSTRATING  THE  USE  OF  THE  LOAD 

CELL  PLATFORM 


FIGURE  4.  ILLUSTRATION  OF  STRAIN  GAGE  LOCATIONS  ON  THE  B-727  N40  LANDING  GEAR 


Xr  ■  Front  and  Back  Axle 
\  -  Front  and  Back  Axle 
It  -  Top  and  Bottom  Axle 
*  Top  and  Bottom  Axle 
Y  •  Front  and  Back  of  Strut 
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FIGURE  6.  TOP  DIAGRAM  OF  METRAPLEX  HYBRID  HIGH-SPEED  DATA  ACQUISITION  SYSTEM 
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FIGURE  7.  EXAMPLE  TO  ILLUSTRATE  EFFECT  OF  THE  FILTER;  Fc  =  30  HZ 
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FIGURE  8.  EXAMPLE  TO  ILLUSTRATE  EFFECT  OF  THE  FILTER;  Fc  =  5  HZ 
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AD  =  Takeoff  Distance 
PR^  Maximum  Pitch  Rate 
Vnax  =  Maximum  Ground  Speed 
0inM=  Maximum  Pitch  Angle 
Nx  =  Maximum  Longitudinal  Accel. 
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FIGURE  9.  ILLUSTRATION  OF  SPECIFIC  TAKEOFF  PARAMETERS 
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Nz  =  Maximum  Normal  Accel. 
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FIGURE  10.  ILLUSTRATION  OF  SPECIFIC  LANDING  PARAMETERS 
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FIGURE  11.  ILLUSTRATION  OF  SPECIFIC  RUNWAY  EXIT  PARAMETERS 
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FIGURE  12.  ILLUSTRATION  OF  SPECIFIC  BRAKING  PARAMETERS 
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Aircraft  S-Turn 
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Aircraft  Minimum-Radius  Turn 
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At  «  arbitrary  10  sec  interval 
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ADL  =  Nose  Gear  Axle 
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FIGURE  14.  ILLUSTRATION  OL  SPECEFIC  MINIMUM-RADIUS  TURN  PARAMETERS 
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Angle  (deg)  Speed  (mph)  Accel  (g’s) 


TST63006 


Data  Rate  =  30  Hz 
Fc  =  5  Hz 


FIGURE  15.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  TYPICAL  TAKEOFF  EVENT 
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Angle  (deg)  Speed  (mph)  Accel  (g’s) 


TST11408 


Data  Rate  =  30  Hz 
Fc  =  5  Hz 


FIGURE  16.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  TYPICAL  TAKEOFF  EVENT 
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Pitch  Pitch 


Time  (sec) 

GRAPH  17(a) 


Time  (sec) 

GRAPH  17(b) 


FIGURE  17.  PITCH  ANGLE  AND  PITCH  RATE  DATA  OF  TWO  TYPICAL  TAKEOFF  EVENTS  (a.) 

TST63006  AND  (b.)  TST11408 


TABLE  1.  TABLE  OF  SPECIFIC  EVENT  PARAMETERS  FOR  THE  TAKEOFF  EVENTS 
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FIGURE  19.  HISTOGRAM  OF  TAKEOFF  GROUND  SPEED  SPECIFIC  EVENT  VALUES 
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FIGURE  22.  HISTOGRAM  OF  TAKEOFF  TIME  TO  LIFTOFF  SPECIFIC  EVENT  VALUES 
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FIGURE  23.  HISTOGRAM  OF  TAKEOFF  MAXIMUM  PITCH  RATE  SPECIFIC  EVENT  VALUES 
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FIGURE  24.  HISTOGRAM  OF  TAKEOFF  MAXIMUM  PITCH  ANGLE  SPECIFIC  EVENT  VALUES 
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FIGURE  25.  CORRELATION  PLOT  OF  TAKEOFF  DISTANCE  AND  TAKEOFF  SPEED 
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FIGURE  27.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  TYPICAL  LANDING  EVENT 
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FIGURE  28.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  TYPICAL  LANDING  EVENT 
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FIGURE  29.  COMPARISON  OF  MEASURED  AND  CALCULATED  VERTICAL  SHEAR  AT  TOUCHDOWN 
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FIGURE  31.  HISTOGRAM  OF  LANDING  TOUCHDOWN  GROUND  SPEED  SPECIFIC  EVENT  VALUES 
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FIGURE  32.  HISTOGRAM  OF  LANDING  ADJUSTED  TOUCHDOWN  SPEED  SPECIFIC  EVENT  VALUES 
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FIGURE  33.  HISTOGRAM  OF  LANDING  NORMAL  ACCELERATION  SPECIFIC  EVENT  VALUES 


HISTOGRAM 

Landing  Data 


TD  V-Shear  (kips) 


FIGURE  34.  HISTOGRAM  OF  LANDING  TOUCHDOWN  VERTICAL  SHEAR  SPECIFIC  EVENT  VALUES 
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FIGURE  35.  CORRELATION  PLOT  OF  PITCH  ANGLE  AND  AIRCRAFT  WEIGHT  AT  TOUCHDOWN 


Peak  Nz  (g’s) 


FIGURE  36.  CORRELATION  PLOT  OF  VERTICAL  SHEAR  AND  NORMAL  ACCELERATION 

AT  TOUCHDOWN 
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FIGURE  37.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  TYPICAL  RUNWAY  EXIT  AT  40  KNOTS 
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FIGURE  38.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  TYPICAL  RUNWAY  EXIT  AT  50  KNOTS 
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FIGURE  39.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  TYPICAL  RUNWAY  EXIT  AT  60  KNOTS 
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FIGURE  40.  COMPARISON  IN  TIME  OF  SIDE  FORCE  COEFFICIENT  AND  LATERAL  ACCELERATION 

FOR  (a.)  TST62904  AND  (b.)  TST62907 
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FIGURE  41.  CORRELATION  IN  TIME  OF  SIDE  SHEAR  AND  LATERAL  ACCELERATION 
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FIGURE  43.  COMPARISON  OF  MEASURED  AND  CALCULATED  AVERAGE  SIDE  SHEAR 
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FIGURE  44.  CORRELATION  PLOT  OF  EXIT  GROUND  SPEED  AND  LATERAL  ACCELERATION 


FIGURE  45.  CORRELATION  PLOT  OF  SIDE  SHEAR  AND  LATERAL  ACCELERATION 
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FIGURE  46.  CORRELATION  PLOT  OF  SIDE  SHEAR  AND  EXIT  GROUND  SPEED 
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FIGURE  47.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  HEAVY  BRAKING  EVENT 
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FIGURE  51.  COMPARISON  OF  MEASURED  AND  CALCULATED  AVERAGE  LONGITUDINAL 

ACCELERATION 
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FIGURE  52.  COMPARISON  OF  MEASURED  AND  CALCULATED  AVERAGE  DRAG  SHEAR 
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FIGURE  53.  CORRELATION  PLOT  OF  AVERAGE  DRAG  SHEAR  AND  LONGITUDINAL 

ACCELERATION 
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FIGURE  57.  COMPARISON  IN  TIME  OF  SIDE  FORCE  COEFFICIENT  AND  LATERAL  ACCELERATION 

FOR  (a.)  TST7103  AND  (b.)  TST63008 
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FIGURE  58.  SIDE  SHEAR  AND  AXLE  DIFFERENTIAL  LOAD  PLOTTED  IN  TIME  FOR  THE  TWO 

MAIN  LANDING  GEAR;  40-KNOT  S-TURN 
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FIGURE  59.  SIDE  SHEAR  AND  AXLE  DIFFERENTIAL  LOAD  PLOTTED  IN  TIME  FOR  THE  TWO 

MAIN  LANDING  GEAR;  60-KNOT  S-TURN 
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FIGURE  60.  SIDE  SHEAR  AND  AXLE  DIFFERENTIAL  LOAD  PLOTTED  IN  TIME  FOR  THE  TWO 

MAIN  LANDING  GEAR;  80-KNOT  S-TURN 
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FIGURE  61.  CORRELATION  IN  TIME  OF  SIDE  SHEAR  AND  LATERAL  ACCELERATION 
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FIGURE  63.  COMPARISON  OF  MEASURED  AND  CALCULATED  AVERAGE  SIDE  SHEAR 
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FIGURE  64.  CORRELATION  PLOT  OF  SIDE  SHEAR  AND  LATERAL  ACCELERATION 
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FIGURE  66.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  TYPICAL  MINIMUM-RADIUS  TURN;  CCW 
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FIGURE  67.  MULTIPLE  TIME  TRACE  DATA  PLOT  OF  A  TYPICAL  MINIMUM-RADIUS  TURN;  CW 
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FIGURE  68.  PHOTO  OF  NOSE  GEAR  DURING  MINIMUM-RADIUS  TURN 


FIGURE  69.  PHOTO  OF  RUBBER  TIRE  TRACKS  DEPOSITED  ON  PAVEMENT 
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FIGURE  70.  COMPARISON  IN  TIME  OF  SIDE  SHEAR  AND  AXLE  DIFFERENTIAL  LOAD  FOR  (a.) 

TEST32  AND  (b.)  TEST34 
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FIGURE  72.  CORRELATION  PLOT  OF  NOSE  GEAR  AVERAGE  SIDE  SHEAR  AND  AXLE 

DIFFERENTIAL  LOAD 
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APPENDIX  B 
Force  Calibration  Procedure 


The  aircraft  was  placed  on  screw  jack  stands,  supported  at  standard  jack  points,  and  a  zero 
force  output  was  recorded  for  the  three  standard  directions;  vertical  shear,  drag  shear,  and  side 
shear  for  each  landing  gear. 

To  Calibrate  Vertical  Shear: 

Individually,  the  nose  gear,  left  main  gear,  and  right  main  gear  of  the  aircraft  were  lowered 
incrementally  onto  a  load  cell  platform.  Load  cell  output  and  strain  gage  output  were 
recorded  for  several  different  loadings,  up  to  full  aircraft  weight. 

To  Calibrate  Drag  Shear: 

Individually,  the  nose  gear,  left  main  gear,  and  right  main  gear  of  the  aircraft  were  lowered 
onto  a  load  cell  platform  to  approximate  normal  vertical  load.  The  load  cell  platform  was 
loaded  longitudinally  at  several  different  weights  up  to  5000  lbs.  Load  cell  output  and  strain 
gage  output  were  recorded  for  the  different  loadings. 

To  Calibrate  Side  Shear: 

Individually,  the  nose  gear,  left  main  gear,  and  right  main  gear  of  the  aircraft  were  lowered 
onto  a  load  cell  platform  to  approximate  normal  vertical  load.  The  load  cell  platform  was 
loaded  laterally  at  several  different  weights  up  to  5000  lbs.  Load  cell  output  and  strain  gage 
output  were  recorded  for  the  different  loadings. 


APPENDIX  C 
Listing  of  Data  Channels 


Channel  Units  Range 


CG  Acceleration  -  Lateral 
CG  Acceleration  -  Longitudinal 
CG  Acceleration  -  Normal 
Throttle  Position  -  #3  Engine 
Flap  Position 
Elevator  Position 
Rudder  Position 
Aileron  Position 

Brake  Pressure  -  LMG;  Left  Axle 
Brake  Pressure  -  LMG;  Right  Axle 
Brake  Pressure  -  RMG;  Right  Axle 
Brake  Pressure  -  RMG;  Left  Axle 
Ground  Distance 

Ground  Speed;  from  Optical  Sensor 
Strain  Gages 

Vertical  Shear  -  RMG;  Right  Axle 
Vertical  Shear  -  RMG;  Left  Axle 
Vertical  Shear  -  LMG;  Right  Axle 
Vertical  Shear  -  LMG;  Left  Axle 
Vertical  Shear  -  NG;  Right  Axle 
Vertical  Shear  -  NG;  Left  Axle 
Drag  Shear  -  RMG;  Right  Axle 
Drag  Shear  -  RMG;  Left  Axle 
Drag  Shear  -  LMG;  Right  Axle 
Drag  Shear  -  LMG;  Left  Axle 
Drag  Shear  -  NG;  Right  Axle 
Drag  Shear  -  NG;  Left  Axle 
Side  Shear  -  RMG 
Side  Shear  -  LMG 
Side  Shear  -  NG 
'Pitch  Angle 
'Roll  Angle 
'True  Heading 


G’s 

-.5  to  .5 

G’s 

-1  to  1 

G’s 

-1  to  1 

%  travel 

0  to  100 

Degrees 

0  to  40 

Degrees 

-10  to  20 

Degrees 

-25  to  25 

Degrees 

-20  to  20 

psi 

0  to  3000 

psi 

0  to  3000 

psi 

0  to  3000  . 

psi 

0  to  3000 

feet 

121.92  pulses/ft 

mph 

0.2  to  250 

pounds 

pounds 

0  to  200000  (sum) 

pounds 

pounds 

0  to  200000  (sum) 

pounds 

pounds 

0  to  20000  (sum) 

pounds 

pounds 

0  to  50000  (sum) 

pounds 

pounds 

0  to  50000  (sum) 

pounds 

pounds 

0  to  10000  (sum) 

pounds 

0  to  50000 

pounds 

0  to  50000 

pounds 

0  to  10000 

degrees 

-90  to  90 

degrees 

-180  to  180 

degrees 

-180  to  180 

Only  Utilized  for  Takeoffs  and  Landings 

**  Aircraft  instrumentation  uses  standard  sign  convention;  right  handed  coordinates  with 
forward  being  positive  x.  Positive  accelerations  give  negative  forces. 
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APPENDIX  D 


Data  Reduction  Software 
Program  Listing 


Data  Reduction  Software  -  Program  Listing 


'************************************* 
/  *  * 

'  *  PROGRAM  TO:  * 
'  *  Read  signals  from  a  Eidel  binary  file  system  * 
'  *  and  write  to  ASCII  file.  * 
/  *  * 

'  *  Bill  Cavage  * 
'  *  Based  on  a  program  by  Gordon  Hayhoe  * 
/  *  *  *  *************************  *  ******** 


'  declare  subs  and  functions 
/ 

DECLARE  SUB  LPTB  (MM% ,  T!,  B!,  Al ! ( ) ,  A2 ! ( ) ,  BZERO ! ) 

DECLARE  SUB  TTRAN  (Al !  ( ) ,  A2  !  (  )  ,  BO  !  ( ) ,  B1 !  ( ) ,  B2  !  ( )  ,  M%(  T!,  FREQ!,  ABZ ! , 
PHS  !  ) 

DECLARE  SUB  SETTANFILT  () 

DECLARE  SUB  TAN2FILT  ( W ! ,  J&) 

DECLARE  SUB  SETFILTERS  () 

DECLARE  FUNCTION  RIO!  (X!) 

DECLARE  FUNCTION  GETDOUBLEHEX !  (ISIG&,  I&) 

DECLARE  FUNCTION  GETMSDLSD!  (J&,  ISIG&,  IREC&,  SCRATE$ ,  PCHRAT ! ,  ROLRAT ! , 
YAWRAT ! ) 

DECLARE  FUNCTION  GETDEC !  (ISIG&,  IREC&) 

.DECLARE  FUNCTION  GETREAL !  (ISIG&,  IREC&) 

DECLARE  FUNCTION  GETHEX  (ISIG&,  IREC&) 


'  define  variables,  dimension  arrays,  and  set  initial  values 
/ 

DEFLNG  I-N 
DEFSTR  S 

DIM  NADDC(IOO) ,  NADDT (100) ,  NELMC(IOO) 

DIM  CPHDG(IOO) ,  RATE (200) ,  FORCECAL ( 2 0 ) ,  ISIGLINE ( 200 ) ,  IWORD(200), 

SUNIT$ (200 ) 

DIM  SWNAME$ (200) ,  WNUM(200),  WTYP(200),  WFNUM(200),  EIDELSIG ( 200 ) ,  NTAG(200) 
DIM  SKIP$ (200) ,  GAIN (200) ,  OSET(200),  CADD(200),  VALUE(200),  STAT$(200) 

DIM  RMAXVAL ( 2 0 0 ) ,  RMINVAL (200) ,  RMAXTIM (200) ,  RMINTIM(200) ,  TOTVAL(200) 
COMMON  SHARED  TRUE,  FALSE 

COMMON  SHARED  TSP,  ISIGLINEO,  TSIGNAME$ ( ) ,  IWORD ( ) ,  LAST$ ( ) ,  INDX() 

TRUE  =  -1:  FALSE  =  0 

NCNT  =  0:  NCORR  =  1:  SWCH$  =  "OFF" 

CLS  :  PRINT  :  PRINT 


'  set  variable  inputs 

i 

SPATH$  =  "C:\DEC315\JUNE2\" 

STESTSIG$  =  "B-727 " 

STESTNAM$  =  "TEST34” 

SFILNAM$  =  "C:\QB\GLDATA\T34" 

' SFILNAM$  =  "CONS:" 

SDATNAM$  =  " GLDATMAN " 

SCALNAM$  =  "GLCALDAT" 

TANF I LTFREQ  =  5 ! 

TSTART  =  0:  TEND  =90 

SCOND$  =  "Minimum  Radius  Turn" 

STITLE1$  =  "B-727  GROUND/FLIGHT  TEST  DATA  FILE" 
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STITLE2$  =  SCOND$  +  "  Data  from  Manning  Data  Set;  " 
NPRNT  =2 
STATOP$  =  "NO" 

SCRATE$  =  "NO" 

SCOLUM$  =  "NORMAL" 

SCOLTIT$  =  " " 


'  set  file  names 
/ 

SIGFILE$  =  SPATH$  + 
SFORFIL$  =  SPATH$  + 
SDATFIL$  =  SPATH$  + 
SPLBFIL$  =  SPATH$  + 
SCALDAT$  =  SCALNAM$ 
SD ATMAN $  =  SDATNAM$ 
SOPFIL$  =  SFILNAM$ 
STATFIL$  =  SFILNAM$ 
STITLE2$  =  STITLE2$ 
PRINT  " 

PRINT  :  PRINT 


STESTSIG$  +  ".SIG" 
STESTSIG$  +  ".FOR" 
STESTNAM$  +  ".DAT" 
STESTNAM$  +  " . PLB " 

+  ".DAT" 

+  ".DAT" 

■  ".DAT" 

+  ".STA" 

+  SDATFIL$ 

CONVERTING  SDATFIL$ 


'  data  manipulation  input  file 
/ 

OPEN  SD ATMAN $  FOR  INPUT  AS  #5 
INPUT  #5,  NCOLELM 
FOR  I  =  1  TO  NCOLELM 
INPUT  #5,  NELMC(I) 

NEXT  I 

INPUT  #5,  SFCHK$ 

IF  SFCHK$  <>  "#"  THEN 

PRINT,  "  ERROR:  BAD  DATA  FILE  SETUP;  CAN  NOT  INTERPRET 

PRINT  "  S D ATMAN $ ;  "  FOR  MAIN  MODULE" 

PRINT  :  PRINT. 

STATOP$  =  "NO" 

GOTO  500 
END  IF 

INPUT  #5,  NCOLADD 
FOR  I  =  1  TO  NCOLADD 

INPUT  #5,  NADDC(I),  NADDT(I) 

NEXT  I 
CLOSE  #5 


'  check  if  input  file  is  valid 
/ 

OPEN  SDATFIL$  FOR  BINARY  AS  #1 

GET  #1,  ,  ID  '  signed  16  bit  integer 

IF  EOF ( 1 )  =  -1  THEN 

PRINT  "  ERROR:  FILE  DOES  NOT  EXIST." 

PRINT  :  PRINT 
CLOSE  #1 
KILL  DATFILE$ 

STATOP$  =  "NO" 

GOTO  500 
END  IF 
CLOSE  #1 


'  obtain  information  from  format  file 
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OPEN  SFORFIL$  FOR  INPUT  AS  #2 
FOR  I  =  1  TO  9 

INPUT  #2,  TEMP 
NEXT  I 

INPUT  #2,  NWORDS 
INPUT  #2,  TEMP 
INPUT  #2,  TEMP 
INPUT  #2,  SYNCVAL$ 

CLOSE  #2 


'  convert  binary  sync  value  to  decimal 
/ 

AMOUNT  =  0:  VALUE  =  0:  DIGIT  =  0 
FOR  I  =  1  TO  12 

SDIG$  =  MID$ ( SYNCVAL$ ,  13  -  I,  1) 
DIGIT  =  VAL ( SDIG$ ) 

IF  DIGIT  =  1  THEN 

AMOUNT  =  AMOUNT  +  2  A  (I  -  1) 

END  IF 
NEXT  I 

FSYNCVAL  =  AMOUNT 


'  obtain  time  information  from  .plb  file 
/ 

OPEN  SPLBFIL$  FOR  INPUT  AS  #2 

LINE  INPUT  #2,  STEMP$ :  LINE  INPUT  #2,  STEMP$ 

INPUT  #2,  DAYOFMONTH $ 

INPUT  #2,  TIMEOFDAY$ 

FOR  I  =  1  TO  8 

INPUT  #2,  WORDLEN 
NEXT  I 

LINE  INPUT  #2,  STEMP$ :  LINE  INPUT  #2,  TEMP$ 

INPUT  #2,  STEMP$ 

LINE  INPUT  #2,  STEMP$ 

INPUT  #2,  BITFREQ 

FFREQ  =  1000!  *  BITFREQ  /  (WORDLEN  *  NWORDS) 

CLOSE  #2 

TINTERVAL  =  1!  /  FFREQ 
TSP  =  TINTERVAL 
FFNYQUIST  =  1 !  /  (TSP  *  2 ! ) 

IF  TANF I LTFREQ  >  FFNYQUIST  THEN  TANFILTFREQ  =  FFNYQUIST 
T VAL ID  =  TSTART  +  (10  *  TINTERVAL) 


'  calculate  filter  weights,  set  filters 

I 

CALL  SETFILTERS 
CALL  SETTANFILT 


'  determine  word  name,  word  number,  frame  number,  and  word  type 
/ 

NWORDS  =  0 
GNUM  =  0 

OPEN  SIGFILE$  FOR  INPUT  AS  #1 
OPEN  SCALDAT$  FOR  INPUT  AS  #4 
LINE  INPUT  #4,  STEMP$ 

IF  STEMP$  <>  "$"  THEN 

PRINT  "  ERROR:  BAD  CALIBRATION  FILE" 

PRINT  :  PRINT 
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STATOP$  =  "NO" 

GOTO  500 
END  IP 

LINE  INPUT  #1,  STEMP$ 

DO 

LINE  INPUT  #1,  ST$ 

IF  ST$  =  "#"  THEN  EXIT  DO 
IST1  =  INSTR ( ST$ ,  "*«)  +  1 

IST2  =  INSTR (ST$,  ";") 

IST3  =  IST2  -  IST1 

SNFIND$  =  MID$(ST$,  ISTCNT  +1,  1) 
IF  SNFIND$  =  THEN 

SFNUM$  =  "0" 

ELSE 


clear  first  line  of  .SIG  file, 
loop  until  no  more  word  entries 
one  complete  line  from  ,SIG  file. 

number  of  characters  to  " 
number  of  characters  to  first  ; 
length  of  channel  name 


IST4  =  INSTR (ISTCNT  +  1,  ST$ , 

SFNUM$  =  MID$ ( ST$ ,  ISTCNT  +  1,  IST4) 

END  IF 

ISTCNT  =  IST2 
ICNT  =  1 
SWCH$  =  "OFF" 

DO 

ISTCNT  =  ISTCNT  +  1 

SNFIND$  =  MID$(ST$,  ISTCNT,  1) 

IF  SNFIND$  =  " ; "  THEN 
ICNT  =  ICNT  +  1 
IF  ICNT  =  2  THEN 

ISTCNT  =  ISTCNT  +  1 
IST7  =  INSTR (ISTCNT,  ST$, 

IF  IST7  -  ISTCNT  =  1  THEN 

SWNUM$  =  MID$(ST$,  ISTCNT,  1) 

ELSEIF  IST7  -  ISTCNT  =  2  THEN 

SWNUM$  =  MID$(ST$,  ISTCNT,  2} 

ISTCNT  =  ISTCNT  +  1 
ELSEIF  IST7  -  ISTCNT  =  3  THEN 

SWNUM$  =  MID$(ST$,  ISTCNT,  3) 

ISTCNT  =  ISTCNT  +  2 
END  IF 
END  IF 

IF  ICNT  =  6  THEN 

ISTCNT  =  ISTCNT  +  1 

SWTYP$  =  MID$(ST$,  ISTCNT,  1) 

END  IF 

IF  ICNT  =  7  THEN 

ISTCNT  =  ISTCNT  +  1 

IST5  =  INSTR (ISTCNT,  ST$, 

SGAN$  =  MID$(ST$,  ISTCNT,  IST5  -  ISTCNT) 

ISTCNT  =  IST5 

ISTCNT  =  ISTCNT  +  1 

IST6  =  INSTR (ISTCNT,  ST$ , 

SOST$  =  MID$(ST$,  ISTCNT,  IST6  -  ISTCNT) 
ISTCNT  =  IST6 
ICNT  =  ICNT  +  1 
END  IF 

IF  ICNT  =11  THEN 

ISTCNT  =  ISTCNT  +  1 

STAT$  =  MID$ ( ST$ ,  ISTCNT,  8) 

SWCH$  =  "ON" 

END  IF 
END  IF 

LOOP  WHILE  SWCH$  =  "OFF" 


'  presently  not 
'  using  this 
'  routine 
'  to  get 
'  calibration 
'  data  from  .sig 
'  file;  it  has 
'  not  been  correct 
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NUMB  =  VAL ( SWNUM$ )  +  1 
IF  NUMB  >  GNUM  THEN 
GNUM  =  NUMB 
END  IF 

SWNAME$ (NUMB)  =  MID$(ST$,  IST1,  IST3 ) 

WNUM (NUMB)  =  VAL ( SWNUM$ ) 

WTYP ( NUMB )  =  VAL ( S WTYP $) 

WFNUM ( NUMB )  =  VAL(SFNUM$) 

SDIG4$  =  MID$(STAT$,  4,  1) 

NT AG ( NUMB )  =  VAL ( SDIG4  $ ) 

INPUT  #4,  GAIN (NUMB) 

INPUT  #4,  OSET (NUMB) 

INPUT  #4,  SUNIT$ (NUMB) 

NWORDS  =  NWORDS  +  1 

LOOP 

CLOSE  #1 

INPUT  #4,  STEMP$ 

IF  STEMP$  <>  "#"  THEN 

PRINT  "  ERROR:  BAD  CALIBRATION  FILE" 

PRINT  :  PRINT 
STATOP$  =  "NO" 

GOTO  500 
END  IF 
CLOSE  #4 


'  routine  to  eliminate  blank  spaces  in  .sig  file 
/ 

1  =  1 
DO 

IF  SWNAME$(I)  =  ""  THEN 
FOR  J  =  I  TO  GNUM 
WNUM(J)  =  J  -  1 
SWNAME  $ ( J )  =  SWNAME$ ( J  +  1) 

WTYP(J)  =  WTYP ( J  +  1) 

WFNUM (J)  =  WFNUM ( J  +  1) 

GAIN ( J )  =  GAIN ( J  +  1) 

OSET ( J)  =  OSET ( J  +  1) 

NTAG(J)  =  NT AG ( J  +  1) 

SUNIT$(J)  =  SUNIT$(J  +  1) 

NEXT  J 
END  IF 

IF  SWNAME$ ( I )  <>  ""  THEN 
1  =  1  +  1 
END  IF 

LOOP  UNTIL  I  >  NWORDS 


'  mark  unwanted  columns;  mark  columns  to  be  combined 
/ 

FOR  I  =  1  TO  NWORDS 

FOR  J  =  1  TO  NCOLELM 

IF  NELMC(J)  =  I  THEN 
SKIP$ ( I )  =  "YES" 

END  IF 
NEXT  J 

FOR  J  =  1  TO  NCOLADD 

IF  NADDC(J)  =  I  THEN 
CADD ( I )  =  J 
CADD ( NADDT ( J ) )  =  J 
END  IF 
NEXT  J 
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FOR  J  =  1  TO  NCOLADD 

IF  NADDT(J)  =  I  THEN 
SKIP$(I)  =  " YES " 
END  IF 
NEXT  J 
NEXT  I 


'  eliminate  channels  not  tagged 

t 

NCHANELM  =  0 
1  =  1 
DO 

IF  NTAG(I)  =  0  THEN 

NCHANELM  =  NCHANELM  +  1 
FOR  J  =  I  TO  NWORDS 
WNUM(J)  =  J  -  1 
SWNAME$ ( J)  =  SWNAME$(J  +  1) 
WTYP { J)  =  WTYP ( J  +  1) 
WFNUM(J)  =  WFNUM ( J  +  1) 
GAIN ( J)  =  GAIN ( J  +  1) 
OSET(J)  =  OSET ( J  +  1) 
NTAG(J)  =  NT AG ( J  +  1) 
SKIP$(J)  =  SKIP$(J  +  1) 
CADD(J)  =  CADD ( J  +  1) 
NTAG(J)  =  NT AG ( J  +  1) 
SUNIT$(J)  =  SUNIT$ ( J  +  1) 
NEXT  J 
END  IF 

IF  NTAG(I)  <>  0  THEN 
1  =  1  +  1 
END  IF 

LOOP  WHILE  SWNAME$ ( I )  <>  " " 


'  routine  to  eliminate  spaces  form  column  names 
/ 

FOR  I  =  3  TO  NWORDS  -  NCHANELM 
NSTRLEN  =  LEN(SWNAME$ (I) ) 

FOR  J  =  1  TO  NSTRLEN 

SLETSTR$  =  MID$ (SWNAME$ (I) ,  J,  1) 

IF  SLETSTR$  =  "  "  THEN 

MID$ (SWNAME$ (I) ,  J,  1)  = 

END  IF 
NEXT  J 
NEXT  I 


'  open  output  file,  label  top,  and  print  out  headings 

DELTAT  =  TINTERVAL  *  NPRNT 
DFFREQ  =  1  /  DELTAT 
OPEN  SOPFIL$  FOR  OUTPUT  AS  #3 
PRINT  #3,  STITLE1$ 

PRINT  #3,  STITLE2$ 

PRINT  #3,  "Time  Start  =  TSTART ;  Time  Finish  =  "/TEND 

PRINT  #3,  "Time  Between  Samples  =  ";  DELTAT;  DFFREQ;  "Hz  Data  Frequency) 

PRINT  #3, 

PRINT  #3,  "TIME  . 

SWCH$  =  "OFF" 

IF  SCOLUM$  =  "SPECIFIED"  THEN 
PRINT  #3,  SCOLTIT$ 
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ELSE 

FOR  I  =  3  TO  NWORDS  -  NCHANELM 
IF  SKIP$ ( I )  <>  "YES"  THEN 
PRINT  #3,  SWNAME$ ( I ) ;  " 
END  IF 
NEXT  I 

IF  SCRATE$  =  "YES"  THEN 

PRINT  #3,  "ROLRAT  PCHRAT " ; 
END  IF 
END  IF 
PRINT  #3, 

PRINT  #3, 


'  open  Eidel  file  and  find  first  data  word 
/ 

OPEN  SDATFIL$  FOR  BINARY  AS  #1 
ISTART  =  1 

FOR  I  =  ISTART  TO  32000000 
W  =  GETDEC ( 1 ,  I) 

IF  W  =  FSYNCVAL  THEN 
ISTART  =  I 
EXIT  FOR 
END  IF 
NEXT  I 

FOR  I  =  ISTART  +  1  TO  5000 
W  =  GETDEC (1,  I) 

IF  W  =  FSYNCVAL  THEN 

NSIGS  =  {(I  -  ISTART)  \  2) 

EXIT  FOR 
END  IF 
NEXT  I 

' IF  NSIGS  <>  NWORDS  -  NCHANELM  THEN 

'  PRINT  #3 ,  "  ALARM:  DISCREPENCY  IN  NUMBER  OF  WORDS  SPECIFIED1* 

'  PRINT  #3,  "  CHECK  FOR  PROPER  .SIG  FILE" 

NWORDS  =  NSIGS 

'  ELSE 

NWORDS  =  NSIGS 
' END  IF 


'  find  proper  frame  to  match  TSTART 
/ 

SWCH$  =  "OFF" 

TIME  =  0 
I  =  ISTART 
DO 

W  =  GETDEC (1,  I) 

IF  W  =  FSYNCVAL  THEN 

TIME  =  TIME  +  TINTERVAL 
END  IF 
1  =  1  +  1 

LOOP  UNTIL  TIME  >  TSTART 
TIME  =  TIME  -  TINTERVAL 
1  =  1-1 
ISIG  =  I 


'  main  data  output  do  loop 
/ 

DO 
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check  for  sync  word;  search  if  necessary;  be  sure  file  has  data 

W  =  GETDEC ( 1 ,  I) 

IF  W  <>  FSYNCVAL  THEN 
1  =  1  +  1 

IF  EOF ( 1 )  =  -1  THEN 
PRINT  #3, 

PRINT  #3,  *  FILE  OUT  OF  DATA" 

PRINT  "  FILE  OUT  OF  DATA" 

GOTO  500 
END  IF 
PRINT  #3 , 

PRINT  #3,  "  ERROR:  LOSS  OF  SYNC;  CAN'T  FIND  SYNC  VALUE" 

PRINT  #3,  "  SEARCHING  FOR  SYNC  VALUE" 

PRINT  #3, 

DO 

W  =  GETDEC (1,  I) 

IF  W  =  FSYNCVAL  THEN 

PRINT  #3,  "  FOUND  SYNC  WORD ;' RESUMING  EXTRACTION" 

PRINT  #3, 

ISIG  =  I 
NCNT  =0 
EXIT  DO 

ELSE 

1  =  1  +  1 
NCNT  =  NCNT  +  1 
END  IF 

IF  NCNT. >  10000  THEN 

PRINT  #3,  "  CAN'T  FIND  SYNC  WORD;  ENDING  PROGRAM" 

STATOP$  =  "NO" 

GOTO  500 
END  IF 

LOOP  WHILE  Wo  FSYNCVAL 
END  IF 


go  thru  frame  byte  by  byte  getting  data  based  on  signal  type 
FOR  J  =  1  TO  2 

IF  WTYP(J)  =  0  THEN  ' 

W  =  GETDEC (J,  I) 

ELSE 

PRINT  #3,  "  ERROR:  FIRST  2  FRAME  WORDS  ARE  NOT  SYNC" 

PRINT  #3,  "  VALUES" 

STATOP$  =  "NO" 

GOTO  500 
END  IF 

ISIG  =  ISIG  +  2 
NEXT  J 

FOR  J  =  3  TO  NWORDS  '  start  loop  after  sync  words 

IF  WTYP(J)  =  0  THEN  '  read  and  convert  file  data 

IF  J  >  2  THEN  '  based  on  the  signal  type. 

W  =  GETMSDLSD ( J ,  ISIG,  I,  SCRATE$ ,  PCHRAT,  ROLRAT,  YAWRAT) 

ELSE 

W  =  GETDEC (J,  I) 

END  IF 

ELSEIF  WTYP(J)  =  2  THEN 
W  =  GETDEC (J,  I) 

ELSEIF  WTYP(J)  =  8  THEN 
W  =  GETREAL ( J,  I)  . 

CALL  TAN2FILT (W,  J) 
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ELSEIF  WTYP(J)  =  1  THEN 
W  =  GETHEX ( J,  I) 


ELSEIF  WTYP(J)  =  4 

THEN 

W  =  GETREAL ( J, 

I)  * 

5! 

ELSEIF  WTYP(J)  =  5 

THEN 

W  =  GETREAL (J, 

I)  * 

10! 

'  CALL  TAN2FILT (W,  J) 

ELSEIF  WTYP(J)  =  6 

THEN 

W  =  GETREAL (J, 

I)  * 

10!  -  5 

CALL  TAN2FILT (W,  J) 

ELSE 

PRINT  #3,  "  ERROR:  UNRECOGNIZABLE  SIGNAL  TYPE” 

STATOP$  =  "NO" 

GOTO  500 
END  IF 

EIDELSIG ( J)  =  W 

VALUE ( J)  =  (EIDELSIG ( J)  *  GAIN(J) )  +  OSET(J) 

TOTVAL(J)  =  TOTVAL(J)  +  VALUE (J) 


find  and  store  max  and  min  value  of  every  channel 

IF  TIME  >  TVALID  THEN 

IF  SWCH$  =  "OFF"  THEN 

RMAXVAL ( J )  =  VALUE (J) 

RMAXTIM(J)  =  TIME 
RMINVAL ( J)  =  VALUE (J) 

RMINTIM(J)  =  TIME 

ELSE 

IF  VALUE (J)  >  RMAXVAL (J)  THEN 
RMAXVAL (J)  =  VALUE (J) 

RMAXTIM ( J)  =  TIME 
END  IF 

IF  VALUE (J)  <  RMINVAL (J)  THEN 
RMINVAL (J)  =  VALUE (J) 

RMINTIM ( J)  =  TIME 
END  IF 
END  IF 
END  IF 

ISIG  =  ISIG  +  2 
NEXT  J 
I  =  ISIG 

IF  TIME  >  TVALID  THEN  SWCH$  =  "ON” 


eliminate  unwanted  time  records 

IF  NPRNT  <>  NCORR  THEN 
NCORR  =  NCORR  +  1 
GOTO  200 
END  IF 


sort  thru  data  combine  or  average  columns 

FOR  J  =  3  TO  NWORDS 

IF  CADD(J)  <>  0  THEN 

FOR  K  =  3  TO  NWORDS 

IF  K  <>  J  AND  CADD(J)  =  CADD(K)  THEN 
VALUE (J)  =  VALUE (J)  +  VALUE (K) 
END  IF 
NEXT  K 
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END  IF 
NEXT  J 


'  print  columns  specified  and  derived  values 
/ 

PRINT  #3,  USING  "###.####  TIME; 

FOR  J  =  3  TO  NWORDS 

IF  SKIP$ ( J)  <>  "YES"  THEN 

PRINT  #3,  USING  "####.####  ;  VALUE  ( J)  ; 

END  IF 
NEXT  J 

IF  SCRATE$  =  "YES"  THEN 

PRINT  #3,  USING  "###.####  ••  ;  ROLRAT;  PCHRAT ; 
END  IF 
PRINT  #3 , 

NCORR  =  1 

200 

TIME  =  TIME  +  TINTERVAL 
LOOP  WHILE  TIME  <  TEND 

500 


calculate  mean  and  varience;  print  out  max  and  min 

/ 

IF  STATOP$  =  "YES"  THEN 

OPEN  STATFIL$  FOR  OUTPUT  AS  #6 

PRINT  #6,  "TIME  CHANNEL  STATISTICAL  BREAKDOWN  -  SDATFIL$ 

PRINT  #6,  "Max,  Min,  Median,  and  Average:  ";  SCOND$ 

PRINT  #6,  "Time  Start  =  ";  T START ;  "  :  Time  Stop  =  TEND 

PRINT  #6, 

FOR  I  =  3  TO  NWORDS 

IF  SKIP$ (I)  <>  "YES"  THEN 

RMEDIAN  =  (RMAXVAL(I)  +  RMINVAL  ( I )  •)  /  2 
AVERAGE  =  TOTVAL ( I )  /  (TIME  *  DFFREQ  *  NPRNT) 

PRINT  #6,  SWNAME$ ( I ) ;  ”  -  (”;  SUNIT$(I);  “)" 

PRINT  #6,  "Maximum  Value  =  RMAXVAL ( I ) ;  "  @  time  =  RMAXTIM ( I ) 

PRINT  #6,  "Minimum  Value  =  RMINVAL (I);  »  @  time  =  RMINTIM ( I ) 

PRINT  #6,  "Median  =  " ;  RMEDIAN;  "  :  Average  =  AVERAGE 

PRINT  #6, 

END  IF 
NEXT  I 
CLOSE  #6 
END  IF 


'  close  files  and  end  program 

CLOSE  #3 
CLOSE  #1 
END 
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Subroutines 


DEFSNG  S 

ISIG  =  signal  number  in  binary  file  records,  starting  at  1. 

IREC  —  byte  position  in  binary  file  for  start  byte  of  the  current  record 

FUNCTION  GETDEC  (ISIG,  IREC)  STATIC 
DIM  ID  AS  INTEGER 

GET  #1,  ISIG  *  2  -  2  +  IREC,  ID 
W  =  ID 

IF  ID  <  0  THEN  W  =  W  +  65536! 

GETDEC  =  W  /  16! 

END  FUNCTION 

ISIG  =  signal  number  in  binary  file  records,  starting  at  1. 

'  IREC  =  byte  position  in  binary  file  for  start  byte  of  the  current  record 

FUNCTION  GETDOUBLEHEX  (ISIG,  I) 

DIM  ID  AS  INTEGER,  ID1  AS  INTEGER 
DEFLNG  I-N 

ITEMP  =  ISIG  *  2  +  IREC 
GET  #1,  ITEMP  -  2,  ID 
ID  =  ID  AND  127 
Wl#  =  ID 
ID  =  ID'  \  4 
ID1  =  ID 

Wl#  =  ID  'Wl#'  /  4# 

GET  #1,  ITEMP,  ID 
ID  =  ID  'AND  4095 
ID  =  ID  '\  4 
W2#  =  ID 

IF  ID  <  0  THEN  W2#  =  W2#  +  65536# 

W2#  =  W2#  \  64 

PRINT  SIGN;  ID1;  ID;  Wl#;  W2# 

W2#  =  (Wl#  /  128!  +  W2#  /  131072!)  *  180! 

GETDOUBLEHEX  =  W2# 

END  FUNCTION 

'  ISIG  =  signal  number  in  binary  file  records,  starting  at  1. 

'  IREC  =  byte  position  in  binary  file  for  start  byte  of  the  current  record. 
FUNCTION  GETHEX  (ISIG,  IREC)  STATIC 

DIM  HBO  AS  STRING  *  1,  HB1  AS  STRING  *  1,  HB2  AS  STRING  *  1,  HB3  AS  STRING  * 
DIM  TEMP2  AS  STRING,  ID  AS  INTEGER 
SHARED  FIRSTHEX$ ,  SECONDHEX  AS  STRING 

ITEMP  =  ISIG  *  2  -  2  +  IREC 

GET  #1,  ITEMP  +  0,  HBO:  GET  #1,  ITEMP  +  1,  HB1 
GET  #1,  ITEMP  +  2,  HB2 :  GET  #1,  ITEMP  +  3,  HB3 

ITEMP  =  ASC(HBO):  FIRSTLO  =  ITEMP 

IF  ITEMP  =  0  THEN  TEMP2  =  ” 00"  ELSE  TEMP 2  =  HEX$ (ITEMP) 

IF  LEN (TEMP2 )  =  1  THEN  TEMP2  =  "0"  +  TEMP 2 


'  Signed  16  bit  integer. 


'  Positive  REAL  in  12  bit  range . 


'  signed  16  bit  integer. 

'  positive  REAL  in  16  bit  range, 
positive  REAL  in  12  bit  range. 


ITEMP  =  ASC(HBl):  FIRSTHI  =  ITEMP 

IF  ITEMP  =  0  THEN  FIRSTHEX$  =  "00"  ELSE  FIRSTHEX$  =  HEX$ (ITEMP) 
IF  LEN ( FIRSTHEX$ )  =  1  THEN  FIRSTHEX$  =  "0”  +  FIRSTHEX$ 


FIRSTHEX$  =  FIRSTHEX$  +  TEMP2 


ITEMP  =  ASC (HB2 ) :  SECONDLO  =  ITEMP 
IF  ITEMP  =  0  THEN  TEMP2  =  "00"  ELSE  TEMP2  = 
IF  LEN (TEMP2 )  =  1  THEN  TEMP2  =  "0"  +  TEMP2 


HEX$ (ITEMP) 


ITEMP  =  ASC (HB3 ) :  SECONDHI  =  ITEMP 

IF  ITEMP  =  0  THEN  SECONDHEX  =  "00“  ELSE  SECONDHEX  -  HPrt ittpmdi 
IF  LEN ( SECONDHEX )  =  1  THEN  SECONDHEX  =" 0"+  SECONDHEX  $  ^ 
SECONDHEX  =  SECONDHEX  +  TEMP2  SECONDHEX 


IW1  -  (FIRSTHI  *  256  +  FIRSTLO)  \  64 
IW1  ='  IW1  AND  127 
W1  =  IW1  'ID 

IW2  =  (SECONDHI  *  256  +  SECONDLO)  \  64 
W2  =  IW2 


GETHEX  =  (W1  /  128!  +  W2  /  131072!)  *  180! 
END  FUNCTION 


'  trec  -  hi?oal  !n  ^inar^  file  records,  starting  at  1. 

D^HsGSSf  "( I S IG ^ ^REC  f  1SCRATE$  f  ^ PCHRAT  SROLRAT ^  SSSS?  SKS’ 

SHARED  NWPHDG,  CPHDG ( ) ,  NSIGS,  SWNAME$ ( ) ,  RATE ( ) ,  TINTERVAL 


get  msd/lsd  values  and  apply  formula 

TEMP$  =  SWNAME$(J) 

ITEMP  =  J  *  2  +  IREC 
GET  #1,  ITEMP  -  2,  ID 
Wl#  =  ID 

IF  ID  <  0  THEN  Wl#  =  Wl#  +  65536# 

Wl#  =  Wl#  /  16# 

GET  #1,  ITEMP,  ID 
W2#  =  ID 

IF  ID  <  0  THEN  W2#  =  W2#  +  65536# 

W2#  •=  W2#  /  16# 

W3#  =  (Wl#  *  256!  +  W2#)  *  .0054931 


Signed  16  bit  integer 

Positive  REAL  in  16  bit  range 
Positive  REAL  in  12  bit  range 


'  apply  proper  formula  for  proper  channel 


IF  TEMP$  =  "ROLL _ (MSD) "  THEN 

IF  W3#  >  180#  AND  W3#  <=  360#  THEN  W3#  =  W3# 
calculate  roll  rate 
IF  SCRATE$  =  "YES"  THEN 

IF  ITEMP  >  2  *  NWPHDG  *  NSIGS  THEN 
R  =  0! 

FOR  K  =  1  TO  NWPHDG 

IOFFSET  =  K  *  NSIGS  *  2 

GET  #1,  ITEMP  -  2  +  IOFFSET,  ID 

Wl#  =  ID 

IF  ID  <  0  THEN  Wl#  =  Wl#  +  65536# 
Wl#  =  Wl#  /  16# 

GET  #1,  ITEMP  +  IOFFSET,  ID 
W2#  '=  ID 

IF  ID  <  0  THEN  W2#  =  W2#  +  65536# 


360# 


'  signed  16  bit  integer 


'  change  to  12  bit  range 
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W2#  =  W2#  /  16# 


WPLUS#  =  (Wl#  *  256!  +  W2#)  *  .0054931 

IF  WPLUS#  >  180#  AND  WPLUS#  <=  360#  THEN  WPLUS#  =  WPLUS#  -  360# 

GET  #1,  ITEMP  -  2  -  I OFFSET,  ID  '  signed  16  bit  integer 

Wl#  =  ID 

IF  ID  <  0  THEN  Wl#  =  Wl#  +  65536# 

Wl#  =  Wl#  /  16#  '  change  to  12  bit  range 

GET  #1,  ITEMP  -  IOFFSET,  ID 
W2#  =  ID 

IF  ID  <  0  THEN  W2#  =  W2#  +  65536# 

W2#  =  W2#  /  16# 

WMINUS#  =  (Wl#  *  256!  +  W2#)  *  .0054931 

IF  WMINUS#  >  180#  AND  WMINUS#  <=  360#  THEN  WMINUS#  =  WMINUS#  - 

360# 

R  =  R  +  (WPLUS#  -  WMINUS#)  *  CPHDG(K) 

NEXT  K 
ROLRAT  =  R 
END  IF 
END  IF 

ELSEIF  TEMP$  =  " PITCH_ (MSD) "  THEN 

IF  W3#  >  180#  AND  W3#  <=  360#  THEN  W3#  =  W3#  -  360# 

W3#  =  -W3# 

'  calculate  pitch  rate 
IF  SCRATE$  =  “YES”  THEN 

IF  ITEMP  >  2  *  NWPHDG  *  NSIGS  THEN 
R  =  0! 

FOR  K  =  1  TO  NWPHDG 

IOFFSET  =  K  *  NSIGS  *  2 

GET  #1,  ITEMP  -  2  +  IOFFSET,  ID  '  signed  16  bit  integer 

Wl#  =  ID 

IF  ID  <  0  THEN  Wl#  =  Wl#  +  65536# 

Wl#  =  Wl#  /  16#  '  change  to  12  bit  range 

GET  #1,  ITEMP  +  IOFFSET,  ID 
W2#  =  ID 

IF  ID  <  0  THEN  W2#  =  W2#  +  65536# 

W2#  =  W2#  /  16# 

WPLUS#  =  (Wl#  *  256!  +  W2#)  *  .0054931 

IF  WPLUS#  >  180#  AND  WPLUS#  <=  360#  THEN  WPLUS#  =  WPLUS#  -  360# 

GET  #1,  ITEMP  -  2  -  IOFFSET,  ID  '  signed  16  bit  integer 

Wl#  =  ID 

IF  ID  <  0  THEN  Wl#  =  Wl#  +  65536# 

Wl#  =  Wl#  /  16#  '  change  to  12  bit  range 

GET  #1,  ITEMP  -  IOFFSET,  ID 
W2#  =  ID 

IF  ID  <  0  THEN  W2#  =  W2#  +  65536# 

'  W2#  =  W2#  /  16# 

WMINUS#  =  (Wl#  *  256!  +  W2#)  *  .0054931’ 

IF  WMINUS#  >  180#  AND  WMINUS#  <=  360#  THEN  WMINUS#  =  WMINUS#  - 

360# 

R  =  R  +  (WPLUS#  -  WMINUS#)  *  CPHDG(K) 

NEXT  K 
PCHRAT  =  -R 
END  IF 
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END  IF 

ELSEIF  TEMP$  =  " P/HDG_ (MSD) “  THEN 
W3#  =  W3#  -  PCORRECT# 

IF  SGN (W3#)  =  -1  THEN  W3#  =  W3#  +  360# 

'  calculate  yaw  rate 
IF  SCRATE$  =  "YES"  THEN 

IF  ITEMP  >  2  *  NWPHDG  *  NSIGS  THEN 
R  =  0! 

FOR  K  =  1  TO  NWPHDG 

IOFFSET  =  K  *  NSIGS  *  2 

GET  #1,  ITEMP  -  2  +  IOFFSET,  ID  '  signed  16  bit  integer 

Wl#  =  ID 

IF  ID  <  0  THEN  Wl#  =  Wl#  +  65536# 

Wl#  =  Wl#  /  16#  '  change  to  12  bit  range 

GET  #1,  ITEMP  +  IOFFSET,  ID 
W2#  =  ID 

IF  ID  <  0  THEN  W2#  =  W2#  +  65536# 

W2#  =  W2#  /  16# 

WPLUS#  =  (Wl#  *  256!  +  W2#)  *  .0054931 

IF  SGN (WPLUS#)  =  -1  THEN  WPLUS#  =  WPLUS#  +  360# 

GET  #1,  ITEMP  -  2  -  IOFFSET,  ID  '  signed  16  bit  integer 

Wl#  =  ID 

IF  ID  <  0  THEN  Wl#  =  Wl#  +  65536# 

Wl#  =  Wl#  /  16#  '  change  to  12  bit  range 

GET  #1,  ITEMP  -  IOFFSET,  ID 
W2#  =  ID 

IF  ID  <  0  THEN  W2#  =  W2#  +  65536# 

W2#  =  W2#  /  16# 

WMINUS#  =  (Wl#  *  256!  +  W2#)  *  .0054931 

IF  SGN (WMINUS#)  =  -1  THEN  WMINUS#  =  WMINUS#  +  360# 

R  =  R  +  (WPLUS#  -  WMINUS#)  *  CPHDG(K) 

NEXT  K 
YAWRAT  =  R 
END  IF 
END  IF 

ELSEIF  LEFT$ (TEMP$,  3)  =  "EPR"  THEN 
W3#  =  W3#  *  .008333 
ELSEIF  TEMP$  =  "OPT/DIST_MSD"  THEN 
W3#  =  Wl#  *  16  +  W2#  /  256 
END  IF 

W  =  W3# 

GETMSDLSD  =  W 
END  FUNCTION 

'  ISIG  =  signal  number  in  binary  file  records,  starting  at  1 . 

'  IREC  =  byte  position  in  binary  file  for  start  byte  of  the  current  record. 
FUNCTION  GETREAL  (ISIG,  IREC)  STATIC 
SHARED  NWPHDG,  CPHDG (),  NSIGS 
DIM  ID  AS  INTEGER 

IPOSN  =  ISIG  *  2  +  IREC 

GET  #1,  ISIG  *  2  -  2  +  IREC,  ID  '  signed  16  bit  integer. 

W  =  ID 

IF  ID  <  0  THEN  W  =  W  +65536!  '  positive  REAL  in  16  bit  range. 

GETREAL  =  W  /  65536!  '  put  in  range  0.0  to  1.0. 
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«* 


GET  #1,  ISIG  *  2  -  2  +  IREC,  ID 
W  =  ID 

IF  ID  <  0  THEN  W  =  W  +  65536' 

W  =  W  /  65536! 

GETREAL  =  W 

END  FUNCTION 


signed  16  bit  integer. 

positive  REAL  in  16  bit  range, 
put  in  range  0.0  to  1.0. 


DEFINT  I-N 

SUB  LPTB  (MM,  T,  B,  A1 ( ) ,  A2 ( ) ,  BZERO) 
M  =  MM 


ANG  =  3.141593  *  B  *  T 
FACT  =  SIN (ANG)  /  COS (ANG) 
Ml  =  M  -  M  /  2 
F  =  1! 

FFN  =  M 


SECTOR  =  3.141593  /  FFN 
WEDGE  =  SECTOR  /  2! 

FOR  I  =  1  TO  Ml 
FFN  =1-1 


ANG  =  FFN  *  SECTOR  +  WEDGE 
AM  =  FACT  *  SIN (ANG) 

BM  =  FACT  *  COS (ANG) 

AMS  =  AM  *  AM 
DEN  =  (1!  +  BM)  A  2  +  AMS 
A1(I)  =  -2!  *  ((1!  -  BM  *  BM) 
A2(I)  =  ((1!  -  BM)  A  2  +  AMS) 
F  =  F  *  (1!  +  A1(I)  +  A2 (I) ) 

NEXT  I 

BZERO  =  F  A  (1!  /Ml) 


-  AMS) 
/  DEN 
/  4! 


STATIC 


/  DEN 


END  SUB 


DEFSNG  I-N 

FUNCTION  RIO  (X)  STATIC 
DEFINT  I-N 

Y  =  X  /  2!:  T  =  IE-08:  E  =  1!:  DE  =  1 ' 
1  =  0 
DO 


1  =  1  +  1 

DE  =  DE  *  Y  /  I :  SDE  =  DE  *  DE :  E 
LOOP  UNTIL  E  *  T  >  SDE 
RIO  =  E 
END  FUNCTION 


E  +  SDE 


T 


DEFLNG  I-N 

SUB  SETFILTERS  STATIC 
SHARED  NWPHDG ,  CPHDG ( ) 


calculate  filter  weights. 


PI  =  3.14159:  TWOPI  = 

NWPHDG  =16 

FB  =  2  *  TSP 

'ALPHA  =  2.783 

ALPHA  =  3.395 

'ALPHA  =  5.653 

AI0  =  RIO (ALPHA) 

FOR  K  =  1  TO  NWPHDG 

AK  =  K 

XIO  =  ALPHA  *  SQR ( 1 ! 


2  !  *  PI 

'  convert  from  Hz 
'  35  dB 
'  40  dB 
'  60  dB 


to  range  0  to  0.5 


-  (AK  *  AK)  /  (NWPHDG  *  NWPHDG) ) 
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CPHDG(K)  =  RIO(XIO)  /  AIO 
OMEGA  =  TWOPI  *  AK  *  FB 

TEMP  =  1!  *  (SIN (OMEGA)  /  AK  -  TWOPI  *  FB  *  COS (OMEGA) )  /  (AK  *  PI) 

CPHDG(K)  =  CPHDG(K)  *  TEMP  /  TSP 
NEXT  K 

1=0:  Y  =  0 ! 

FOR  K  =  1  TO  NWPHDG 

Y  -  Y  +  TSP  *  ((I  +  K)  -  (I  -  K) )  *  CPHDG(K)  'Linear  slope  =  TSP  per  step 

1 

NEXT  K 

YSLOPE  =  Y 

FOR  I  =  1  TO  NWPHDG 

CPHDG(I)  =  CPHDG(I)  /  YSLOPE 
NEXT  I 

GOTO  NOCHECK 
FOR  I  =  1  TO  -NWPHDG 
PRINT  I;  CPHDG(I) 

NEXT  I 

FOR  IFR  =  0  TO  18 
FREQ  =  IFR  *  .02 
HF  =  0! 

FOR  I  =  1  TO  NWPHDG 
AI  =  I 

HF  =  2 1  *  CPHDG(I)  *  SIN (TWOPI  *  AI  *  FREQ)  +  HF 

NEXT  I 

PRINT  USING  "###.####  " ;  FREQ  /  TSP;  HF  /  TWOPI 
NEXT  IFR 

PRINT  TSP 
FOR  I  =  0  TO  1 

Y  =  0! 

FOR  K  =  1  TO  NWPHDG 

Y  =  Y  +  TSP  *  ((I  +  K)  -  (I  -  K) )  *  CPHDG(K) 'Linear  slope  =  TSP  per  step 

1 

NEXT  K 

PRINT  I;  Y;  YSLOPE 
NEXT  I 
END 

NOCHECK : 

END  SUB 

SUB  SETTANFILT  STATIC 

DIM  AL1  ( 100  )  ,  AL2  ( 100  )  ,  BLO(IOO),  BLl(100),  BL2(100),  HF(100) 

DIM  MM  AS  INTEGER 

SHARED  B01 ,  B02,  All,  A12 ,  A21,  A22,  TANFILTFREQ 

TWOPI  =  6.2831853#:  PI  =  TWOPI  /  2! 

MM  =  4 

FB  =  TANFILTFREQ 
'  PRINT  TSP 

'  High  pass  section. 

CALL  LPTB (MM,  TSP,  FB,  ALl ( ) ,  AL2 ( ) ,  BZERO) 

FOR  I  =  1  TO  MM  /  2 

BL0(I)  =  1!  *  BZERO:  BL1(I)  =  2  !  *  BZERO:  BL2(I)  =  1!  *  BZERO 

'  PRINT  USING  "###.#####  ALl (I);  AL2(I);  BL0(I);  BL1(I);  BL2(I) 

NEXT  I 
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B01  =  BLO(l):  B02  =  BL0(2) 

All  =  AL1(1):  A12  =  AL1(2):  A21  =  AL2 ( 1 ) :  A22  =  AL2(2) 

GOTO  SKIPCHECK 
FOR  I  =  1  TO  15 
FREQ  =1*2 

CALL  TTRAN (ALIO  ,  AL2 ( ) ,  BLO ( ) ,  BL1 () ,  BL2 ( ) ,  MM,  TSP,  FREQ,  ABZ,  PHS) 
PRINT  USING  "##.##  ###.#####  FREQ;  ABZ 

HF ( I )  =  ABZ 
NEXT  I 
END 

SKIPCHECK: 

END  SUB 

SUB  TAN2FILT  (W,  J)  STATIC 

DIM  WM1 ( 80 ) ,  WM2 (80),  Y1M1(80),  Y1M2(80),  Y2M1(80),  Y2M2(80) 

SHARED  B01 ,  B02,  All,  A12,  A21,  A22 

Y1  =  B01  *  (W  +  2!  *  WM1(J)  +  WM2 ( J) )  -  All  *  Y1M1(J)  -  A21  *  Y1M2 ( J) 

Y2  =  B02  *  (Y1  +  2!  *  Y1M1(J)  +  Y1M2 (J) )  -  A12  *  Y2Ml(J)  -  A22  *  Y2M2(J) 
Y1M2 (J)  =  Y1M1(J):  Y1M1 ( J)  =  Y1 
Y2M2 ( J)  =  Y2M1 ( J) :  Y2M1 ( J)  =  Y2 
WM2 ( J)  =  WM1(J):  WM1(J)  =  W 
W  =  Y2 

END  SUB  ■ 

DEFINT  I-N 

SUB  TTRAN  (A1  (  )  ,  A2  (  )  ,  BOO,  B1  ()  ,  B2  ( )  ,  M,  T,  FREQ,  ABZ,  PHS)  STATIC 

FACT  =  6.283185  *  T 
IP  =  M  -  M  /  2 
ADD  =  0! 

PREV  =  0! 

FD  =  FREQ  *  FACT 

51  =  SIN (FD) 

Cl  =  COS (FD) 

A  =  2 !  *  FD 

52  =  SIN (A) 

C2  =  COS (A) 

ABSA  =1! 

PHSA  =  0! 

FOR  J  =  1  TO  IP 

AR  =  B0 ( J)  +  B1 ( J)  *  Cl  +  B2 (J)  *  C2 
AI  =  -B1 ( J)  *  SI  -  B2 (J)  *  S2 
ANM  =  AR  *  AR  +  AI  *  AI 
PND  =  0 ! 

IF  AI  <>  0!  OR  AR  <>  0!  THEN  PND  =  ATN (AI  /  AR) 

AR  =  1!  +  A1(J)  *  Cl  +  A2 ( J)  *  C2 

AI  =  -A1(J)  *  SI  -  A2(J)  *  S2 

ABSA  =  ABSA  *  ANM  /  (AR  *  AR  +  AI  *  AI) 

NEXT  J 

ABZ  =  SQR(ABSA) 

END  SUB 
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APPENDIX  E 

Axle  Differential  Load  Analysis 


Figure  El  is  a  diagram  of  a  landing  gear 
under  loading  due  to  lateral  motion. 
Intuitively,  this  motion  causing  a  net 
side  force  on  the  landing  gear  will  also 
cause  a  moment  about  the  landing  gear 
axle,  causing  axle  differential  load. 

Axle  differential  load  is  defined  by  the 
following  as  illustrated  in  figure  El. 

ADL  =  FVSr  -  Fk,  (El) 

In  this  equation  ADL  is  axle  differential 
load. 


vs 


y 


From  the  load  diagram  given  as  figure 
E2  it  is  obvious  that  the  total  side  shear 
on  the  landing  gear  is  equal  to  the  sum 
of  the  side  force  reaction  forces  as  given 
below. 


F  = 
rss 


SSI 


+  F, 


SSr 


(E2) 


Side  force  is  manifested  by  lateral 
acceleration  Ny  as  related  by  the 
following  equation. 

Fss  -  N,  WAC  (E3) 
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FIGURE  E2:  LOAD  DIAGRAM  OF  LANDING  GEAR 


By  summation  of  moments  about  the  labeled  moment  center  the  following  equation  can  be 
obtained. 


0  = 


FysftwB 


+  ^SSr  rt  +  ^SSl  rt 
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In  the  above  equation  r,  is  the  loaded  radius  of  the  tire  and  D^  is  the  wheel  base  distance  or 
contact  patch  distance.  Collecting  side  forces  and  vertical  forces  the  following  manipulations 


E-l 


can  be  performed. 


F VsPwB  FVsP 
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F SSr  rt  +  F SSI  rt 


(E5) 
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(Fvs:  FVSr  )  rt(FSSr  +  FSSl  ) 


(E6) 


We  can  substitute  total  side  shear  and  axle  differential  load  as  defined  above  and  solve  for 
Fadl  t0  obtain  the  following  equation. 


ADI 


2rr  FSS 
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For  the  B-727  rt  is  21.5  inches  and  DWB  =  36  inches.  This  gives  the  following  relation. 


ADL 


=  1.2  F, 


ss 


(E8) 
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Measurements  show  that: 

fadl  *  Fss  m 
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